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ABSTRACT 


SYNTHESIS  METHODS  APPLICABLE  TO  SEM 


POLE-PAIR  ADMITTANCES 


The  applicability  of  the  first  layer  conjugate  pole  pair  at  a  time  equiv¬ 
alent  circuit  synthesis  for  energy  collecting  structures  is  considered.  It 
is  found  that,  at  least  for  some  highly  resonant  structures,  the  whole  eigen- 
set  of  poles  can  be  approximately  represented  by  the  first  layer  dominant 
pole.  The  driving  point  admittance  function  can  then  be  synthesized  in  the 
form  of  a  parallel  connection  of  simple  pole-pair  modules.  The  biquadratic 
functions  corresponding  to  each  conjugate  pole  pair  are,  in  general,  not  pos¬ 
itive  real  and  thus  physically  realizable.  It  is  found  that  after  negligible 


resistive  padding  approximate  realization  is  possible.  The  padded  function 

is  a  minimum  biquadratic  -  the  real  part  of  it  is  zero  at  some  frequency  . 

If  hi  =  0  or  w  =  »,  the  function  can  be  easily  realized  by  continued  fraction 
o  o 

expansion.  If  however  uiq  is  finite,  it  is  shown,  that  the  only  transformer¬ 
less  configuration  derivable  is  the  Bott-Duffin  network  .^The  explicit  form 
of  this  network,  as  well  as  Brune  and  Drjalington  network,  is  derived.  It  is 
demonstrated  that  the  Darlington  network  reduces  to  Brune  network  or,  if 
surplus  factors  are  used,  to  Bott-Duffin  network.  The  Miyata  procedure  is 
not  applicable  to  minimum  biquadratics. 

^After  further  approximation  the  Bott-Duffin  network  can  be  reduced  to 
simple  four  element  structure  which  has  been  shown  to  give  favorable  results 
for  the  thin  wire  dipole  and  loop  antennas. 

It  is  demonstrated  that  the  equivalent  networks  derived  are  very  sensi¬ 
tive  to  changes  in  element  values.  This  and  the  influence  of  the  parasitic 
effects  associated  with  real  circuit  elements  can  unfavorably  affect  the 
practicability  of  the  derived  networks,  should  they  be  constructed.'  - 
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Section  I 


INTRODUCTION 


The  need  for  finite,  lumped  equivalent  networks  valid  over  a  significant 
bandwidth  in  the  analysis  and  design  of  antennas  and  scatterers  has  been 
widely  recognized  for  a  long  time.  Such  networks  may  be  useful,  for  example, 
in  experiments  designed  to  test  the  terminal  circuits  intended  to  operate 
with  an  antenna,  before  the  antenna  is  actually  constructed  [35].  Also, 
existing  circuit  analysis  programs  can  be  used  to  calculate  the  frequency  and/ 
or  time  response  of  the  structure  under  consideration. 

The  Singularity  Expansion  Method  (SEM)  is  a  powerful  and  convenient 
means  of  constructing  equivalent  networks  for  antennas  and  scatterers  [4,5, 
29,30,39].  From  the  SEM  data  the  input  admittance1  of  a  single-port  struc¬ 
ture  can  be  expressed  as 


Y(s)  =  l  Y  (s) 
n  n 

where  Y  ,  the  terminal  eigenadmittances ,  are  given  by 
n 


Y 

n 


(s) 


+  Yent(s)  +  Ybl(s) 
n  n 


(i) 


(2) 


where  n  indexes  over  the  eigenmode  sets  of  poles,  the  index  i  distinguishes 
among  the  poles  associated  with  the  given  eigenmode,  and  a^  are  the  resi¬ 
dues  associated  with  the  poles  s  . .  The  yent  and  the  Ybl  are ,  respectively , 

m  n  n 


Admittance  function  resulting  from  the  short-circuit  boundary  value 
problem  is  used  in  this  work.  However,  analogous  results  are  valid  for 
the  impedance  function  derived  from  the  open-circuit  boundary  value 
problem. 


the  possible  entire  function  and  branch  integral  function  contributions. 
Often  it  is  more  convenient  to  rewrite  (2)  in  the  form 


(3) 


Y  (s) 
n 


+  v'ent(s)  +  Y  *bi  (s) 

sni  /  n 


where  the  entire  function  has  been  appropriately  modified  to  compensate  for 
the  constant  terms  included  in  the  summation.  The  terms  under  the  summation 
signs  in  (2)  and  (3)  are  known,  respectively,  as  pole  admittances  and  modi¬ 
fied  pole  admittances  [5]. 

It  should  be  pointed  out  that  for  numerically  or  experimentally  derived 
SEM  descriptions  the  entire  function  and  the  branch  integral  contributions 
are  not  explicitly  identifiable.  To  date,  it  also  appears  that  for  practi¬ 
cal  circuit  synthesis  it  is  necessary  that  both  of  these  terms  be  zero. 
Further,  for  many  problems  of  interest,  they  seem  not  to  occur.  Therefore, 
the  branch  integral  and  entire  function  terms  are  dropped  henceforth  from 
consideration . 

The  SEM  data  (the  poles  and  residues)  for  a  given  structure  can  be 
extracted  in  three  ways:  1)  from  measurements;  2)  from  an  integral  equa¬ 
tion  formulation  by  method  of  moments  technique;  or  for  certain  simple 
structures,  3)  from  an  integral  equation  formulation-  analytically.  As  a 
typical  example,  the  location  of  poles  for  dipole  and  loop  wire  antennas  is 
shown  in  Figures  1  and  2  [42,43,45,7],  The  first  layer  poles  (poles  closest 
to  the  j-axis)  with  the  associated  residues  are  listed  for  both  structures 
in  Table  l1  and  Table  2. 


The  first  pole  of  the  dipole  has  been  slightly  modified  (within  a 
5%  accuracy  margin)  by  Streable  and  Pearson  in  order  to  make  the  modified 
pole-pair  admittance  physically  realizable. 
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TABLE  1 

First  layer  poles  and  residues  for  center-fed  cylindrical 
antenna  with  L/d=100 

L  is  the  length  of  the  dipole  and  d  is  the  diameter  of  the 
wire.  The  even  poles  with  zero  residues  are  omitted. 


No. 

Poles1 

Residues1 

1. 

-.08427  +  j  0.91582 

.1112  x  10-2  +  j  .3121  x  10~3 

3. 

-.1473  +  j  2.870 

.1319  x  10-2  +  j  .3301  x  10-3 

5. 

-.1877  +  j  4.834 

.1423  x  10-2  +  j  .3521  x  10~3 

7. 

-.2177  +  j  6.792 

.1496  x  10"2  +  j  .3699  x  10_3 

9. 

-.2426  +  j  8.736 

.1557  x  10"2  +  j  .3839  x  10~3 

Normalized  to  ctt/L,  c  -  speed  of  light. 

2 

Slightly  modified  by  Streable  and  Pearson. 

The  admittance  (1)  can  be  realized  only  approximately  by  a  finite 
lumped  network  because  the  number  of  modes  is  infinite  and  must  be  truncated. 
Also,  the  number  of  poles  associated  with  each  eigenmode  is  infinite  for 
some  structures  and  only  a  finite  number  of  them  can  be  used.  The  truncated 
admittance  function  may  be  no  longer  positive  real  (PR)1  and  thus  realizable. 

In  such  case  still  useful  approximate  realization  can  often  be  obtained 
by  finding  pole  groupings  which  give  PR  (or  nearly  PR)  admittance  function 
and/or  by  resistive  padding.  Both  approaches  were  first  investigated, 
successfully,  by  Streable  and  Pearson  [39]. 


See  Section  II  for  explicit  definition  of  positive  real  functions. 


n  =  2ln(2ir5) 


TABLE  2 


First  layer  poles  and  residues  of  loop  wire  antenna  with  ft=15. 

ft=21n(2irb/a) ,  where  a  and  b  are  the  radius  of  the  wire  and  the 
radius  of  the  loop,  respectively. 


No. 

Poles^ 

Residues3- 

B 

-.7487  x  10_1  +  j  1.0388 

.5301  +  j  .8934  x  10~4 

2. 

-.1083  +  j  2.0526 

.5726  +  j  .9395  x  10~4 

3. 

-.1340  +  j  3.0625 

.6033  +  j  .9873  x  10_4 

4. 

-.1558  +  j  4.0706 

.6279  +  j  .1033  x  10_3 

5. 

-.1751  +  j  5.0777 

.6489  +  j  .1076  x  10_3 

6. 

-.1927  +  j  6.0840 

.6672  +  j  .1117  x  10~3 

7. 

-.2091  +  j  7.0897 

.6837  +  j  .1157  x  10"3 

8. 

-.2245  +  j  8.0950 

.6988  +  j  .1195  x  10-3 

-.2391  +  j  9.0999 

.7126  +  j  .1231  x  10-3 

-.2529  +  jl0.1046 

.7255  +  j  .1267  x  10~3 

Normalized  to  c/b.  There  is  also  a  pole  at  the  origin  with  the 
real  residue  a  =4.6199  x  10-4.  This  pole  can  be  easily  real¬ 
ized  as  a  single  shunt  inductance  L=l/a„„  and  is  dropped  from 
.  ,  00 
consideration  here. 

One  pole  grouping  possible,  which  leads  to  PR  function  (although 
recent  investigations  indicate  that  the  branch  integral  contribution  must 
be  retained  in  the  terminal  eigenadmittance  for  some  structures)  is  the 
grouping  according  to  eigenmodes  [29].  This,  however,  has  the  disadvantage 
that  the  number  of  poles  increases  dramatically  for  higher  modes  (see 
Figures  1  and  2)  resulting  in  impractically  complex  network.  Also,  the 
third-  and  higher-layer  poles  are  usually  hard  to  extract,  and  for  many 
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even  simple  structures  have  not,  so  far,  been  extracted. 

Even  if  the  original  function  (1)  is  PR,  a  decomposition  of  it  into 
simpler,  realizable  terms  is  desirable.  Such  decomposition  can  simplify  the 
synthesis  procedure  significantly  and  also  can  lead  to  simpler  network  with 
more  desirable  topology. 

One  particularly  attractive  way  of  decomposition  is  breaking  the  origi¬ 
nal  admittance  function  (1)  into  conjugate-pole-pair  modules,  as  first 
suggested  by  Baum  [5]  and  investigated  further  by  Streable  and  Pearson  [39]. 

Because  the  first  layer  poles  are  so  close  to  the  j-axis  (see  Figures 
1  and  2) ,  they  give  the  main  contribution  to  the  input  admittance  of  the 
structure.  The  question  immediately  arises  if  their  contribution  alone 
could  reasonably  good  approximate  the  whole  admittance  over  the  frequency 
range  of  interest.  That  this  is  the  case  (at  least  for  some  high-Q  struc¬ 
tures)  has  also  been  demonstrated  by  Streable  and  Pearson  [39],  In  this 
procedure  the  whole  eigenmode  pole  group  is  represented  by  a  single,  first- 
layer  pole. 

The  first-layer  pole  grouping  is  particularly  useful  in  cases  when  the 
SEM  data  is  determined  from  measurements,  when  the  higher-order  layers  of 
poles  are  extremely  difficult  to  find  [34]. 

The  main  concern  of  this  study  is  the  synthesis  of  the  admittance 
function  (1)  on  the  conjugate-pole-pair  basis.  Thus,  the  realization  of  (1) 
in  the  form 

N 

Y(s)  =  Z  Y. (s)  (4) 

i 


where 


Yi(s) 


(5) 


8 


is  investigated.  In  (5)  the  asterisk  denotes  the  complex  conjugate,  and 
is  the  padding  necessary  to  render  PR.  The  index  i  distinguishes  among 
the  conjugate  pole  pairs. 

After  putting  the  terms  over  common  denominator,  the  admittance  (5) 
takes  the  form  of  a  general  biquadratic  function 


f  (s) 


2 

As  +  BS  +  C 
_  2 

Ds  +  Es  +  F 


(6) 


This  function  received  a  great  deal  of  attention  in  the  literature,  because 
explicit  networks  can  be  derived  for  it.  The  results  scattered  in  the  lit¬ 
erature  [13,14,28,33,20,3,17,41,48,19,36,46,18]  are  used  heavily  in  this 
report.  The  emphasis  of  this  work,  however,  is  on  the  approximate  synthesis 
stemming  from  initially  non-PR  biquadratics  derived  from  SEM  data.  This  has 
important  consequences  on  what  synthesis  procedures  are  applicable. 

If  not  otherwise  stated  it  will  be  assumed  that  f(s)  represents  an 
admittance  function  Y(s). 

The  plan  of  the  work  is  as  follows.  In  Section  II  are  summarized  the 
basic  concepts  concerning  the  PR  functions,  to  which  references  will  be  made 
in  later  parts  of  this  study.  In  Section  III  the  PR-ness  conditions  for 
biquadratic  functions  are  derived,  the  necessary  padding  (if  the  function 
is  non-PR)  is  computed,  and  synthesis  methods  for  some  special  cases  are 
considered. 

In  Section  IV  a  study  on  the  applicability  of  the  various  general 
synthesis  procedures  to  a  minimum  biquadratic  function  is  undertaken. 
Explicit  Brune,  Bott-Duffin  and  Darlington  networks  are  derived.  In  Section 
V  some  important  problems  connected  with  the  practical  realization  of  the 
networks,  such  as  the  influence  of  parasitic  elements,  are  briefly 


considered.  A  computer  program  for  the  synthesis  of  conjugate-pole-pair 
modules,  implementing  the  results  of  previous  sections,  is  included  in  the 
Appendix. 

Throughout  this  study  some  concepts  are  illustrated  with  the  particular 
examples  of  the  straight  wire  and  loop  wire  structures  with  the  parameters 
given  in  the  Figures  1  and  2,  and  in  Tables  1  and  2.  The  particular  values 
L=100m  and  b=100m  are  used  in  synthesis  in  order  to  fix  explicit  values 
for  circuit  elements. 


Section  II 


POSITIVE  REAL  FUNCTIONS  AND  BASIC  CONCEPTS 

A  positive  real  (PR)  function  f(s)  is  defined  as  a  function  of  complex 
variable  s  =  0  +  ju>  satisfying  the  two  requirements 

Re  f  (s)  >_  0  for  a  0, 

,  (7) 

and 

f(s)  real  when  s  is  real. 

The  positive  real  concept  was  introduced  by  Brune  [10],  who  proved  that 
every  driving-point  function  of  a  physical  network  is  PR  and  that  every 
rational  function  that  is  PR  can  be  realized  using  lumped  RLCM  elements.  It 
was  later  proved  by  Bott  and  Duffin  that  any  rational  PR  function  can  be 
realized  as  a  driving  point  immitance  with  RLC  elements  only  [8] . 

Since  it  is  extremely  difficult  to  apply  the  conditions  (7)  directly, 
an  equivalent  set  of  requirements,  which  are  easier  to  check,  is  usually 
used  [15].  Thus,  a  rational  function  f (s)  is  PR  if  the  following  necessary 
and  sufficient  conditions  are  satisfied  (listed  and  carried  out  in  their 
order  of  difficulty  in  testing  them) : 

A.  It  has  no  poles  or  zeros  in  the  right  half-plane 

B.  Any  j-axis  poles  are  simple  and  have  positive  real  residues 

C.  Re  f  ( jw)  5^  0  for  0  <_  co  <_  «  . 

In  an  attempt  to  realize  a  driving-point  immittance  f(s)  it  is  usually 
worthwhile,  before  applying  a  general  synthesis  procedure,  to  check  if 
there  are  any  poles  or  zeros  on  the  j-axis.  Such  poles  and  zeros  can  be 
easily  removed  from  Y(s)  (susceptance  reduction)  and/or  from  Z(s)  =  1/Y(s) 
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(reactance  reduction)  without  destroying  the  Pr-ness  of  the  original  func¬ 
tion.  When  all  the  j-axis  poles  are  removed  from  an  admittance  function 
Y (s) ,  it  is  called  a  minimum-susceptance  function  [1],  In  a  similar  manner 
we  define  a  minimum-reactance  function  as  an  impedance  Z(s)  that  has  no  poles 
on  the  j-axis.  When,  after  applying  these  steps,  a  minimum-susceptance- 
minimum-reactance  function  is  obtained,  still  another  simplification,  the 
so  called  resistance  reduction  is  often  possible.  Namely,  if  the  condition 
C.  above  is  satisfied  with  the  inequality  sign,  a  conductance  G  can  be 
removed  from  Y(s)  (or  resistance  R  from  Z(s))  still  leaving  the  remaining 
function  PR.  The  real  part  of  this  remainder  function  is  zero  at  some  fre¬ 
quency  u)Q,  therefore  it  is  called  a  minimum-real-part  function.  If  the  zero 

of  the  real  part  occurs  at  oj  =  0  or  uA  =  ®,  it  is  also  a  zero  of  the  whole 

o  0 

function  because  the  imaginary  part  of  a  minimum- reactive-minimum- susceptive 
function  is  already  zero  at  these  extreme  frequencies.1  This  zero  of  Y(s) 
(Z(s))  can  be  subsequently  removed  as  a  pole  of  Z(s)  (Y(s))  and  the  whole 
process  of  real-  and  imagninary-part  reduction  can  be  attempted  again  until 
the  so-called  minimum  function  (minimum-reactance,  minimum-susceptance  and 
minimum-real-part)  is  obtained.  In  some  cases  this  Foster  preamble  tech¬ 
nique  can  lead  to  a  complete,  successful  realization  [38].  In  general, 

^Consider  Y(s),  for  example.  Im  Y(jui)  is  an  odd  rational  function  of  w  and 
can  take  one  of  the  forms 

Im  {Y ( jui)  }  =  id  — or  Xm{ Y ( ju>)  >  =  N^U12^  . 

D(0)  )  D  ((*»'  ) 

Assume  that  Y(s)  is  purely  imaginary  at  u  =  0.  Then  if  the  first  form 
above  holds,  Z(s)  will  have  a  removable  pole  at  s  =  0.  If  the  second  form 
holds,  then  Y(s)  has  a  pole  at  s  =  0  which  can  be  removed.  Similar  argu¬ 
ments  hold  for  the  case  when  the  zero  of  the  real  part  occurs  at  infinity. 
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however,  we  end  up  with  a  minimum  function  with  zero  real  part  at  some  finite 
frequency  u  ,  and  one  of  the  general  synthesis  procedures  must  be  employed. 
The  procedures  in  question  are  those  of  Brune  [10],  Darlington  [12],  Bott 
and  Duffin  [8],  and  Miyata  [26]. 


Section  III 


SYNTHESIS  OF  BIQUADRATIC  FUNCTIONS  -  PRELIMINARY 


In  this  and  the  following  sections  we  will  be  mainly  concerned  with  the 
synthesis  of  the  biquadratic  function  (6)  which  can  represent  both  Z(s)  or 
Y (s) ,  as  the  case  might  be. 

The  PR-ness  condition  for  the  immitance  function  (6)  can  be  easily 
derived  [23,1,21].  If  the  coefficients  of  (6)  are  positive,  the  conditions 
A.  and  B.  of  Section  II  are  automatically  satisfied.  Thus,  to  guarantee 
that  f(s)  be  PR,  we  require  that 


Re{f  ( ja>)  }  = 


mlm2  ‘  nin2 

2  2 
m2  -  n2 


s  =  }<i> 


4  2 

ADio  +  (BE -AF -CD)  ai  +  CF 

_  2  2  2  2 
(F  —  Da)  )  +  E  a) 


(8) 


be  nonnegative  for  0<w<_°°  (condition  C.),  where  m  ,  n^  denote,  respectively, 
the  even  and  odd  parts  of  the  numerator,  and  m2,  n2  denote  the  even  and 
odd  parts  of  the  denominator  of  (6) .  This  leads  to  the  following  necessary 
and  sufficient  condition  for  f(s)  to  be  PR: 


(/AF  -  /CD) 1  <  BE. 


(9) 


If  (9)  is  fulfilled  with  the  equality  sign,  f(s)  is  a  minimum  function. 

For  the  important  special  case  where  f(s)  ->-0  at  infinity,  i.e.  A  =  0, 
we  have 

Bs  +  C 


f(s)  = 


(10) 


Ds  +  Es  +  F 


and  the  PR-ness  condition  (9)  now  reads 


BE  -  CD  >  0. 


(11) 
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When  f (o)  =  0,  i.e.  C  =  0,  we  obtain 


f(s)  = 


As  +  Bs 


with  the  condition 


Ds  +  Es  +  F 


BE  -  AF  >  0. 


The  above  conditions  can  be  easily  expressed  in  terms  of  the  poles  and 
residues.  Thus,  for  the  conjugate  pole-pair  admittance 

a  a* 

Y  (s)  =  r  (14) 

n  s  -  S  s  -  s 

n  n 

with  a  =  a  +  and  s  =  -a  +  ju>  we  get  the  special  case  (10)  where 
n  n  n  n  n  n 


B  =  2a 


C  =  2(aa  -  B  w  ) 
n  n  n  n 


D  =  1 


E  *  2a 


F  =  s 


The  PR-ness  condition  (11)  and  the  requirement  that  the  above  coefficients 
be  positive  give  the  following  restriction  on  the  residues 


T3j  i  a 


This  condition  can  be  given  a  simple  geometrical  interpretation,  as  dis¬ 
cussed  by  Guillemin  [15]. 

The  other  special  case  (12)  results  from  the  modified  conjugate-pole- 


pair  module 


a  a  a  a 

-  ,  ,  n  n  .  n  ,  n 

n  s-s  s  s-s*  s 

n  n  n  n 


This  function  approaches  zero  for  s  -*■  0  and  is  particularly  well  suited 


for  the  modelling  of  the  input  admittance  of  simply  connected  objects. 


After  combination  of  terms  in  (16)  we  get 


A  =  2(f3ui  -ao) 
n  n  n  n 


B  =  2a 

n 


2 

a  )  +  48  a  u) 
n  n  n  n 


D  = 


E 


2a 

n 


The  requirement  that  A  and  B  be  positive  leads  to  the  conditions 


6  u  -  a  a  >  0 
n  n  n  n  — 


and 

2  2 

a  (a)  -  a  )  +  28  a  u  >  0. 

n  n  n  n  n  n 


The  PR-ness  condition  (14)  gives 


2  2  2  2 
-a  a  (a  -  3w  )  +  g  tu  (3a  -  u>  )  >  0. 

nnn  n  nnn  n  — 


(17) 


(18) 


(19) 


It  is  interesting  to  point  out  that  the  PR-ness  tests  (15)  and  (17) 
are  mutual ly-exclusivef  which  means  that  a  conjugate-pole-pair  module  can 
be  realized  only  in  one  of  the  two  forms  (modified  or  unmodified) ,  but  not 
in  both. 

As  an  example,  the  test  values  for  the  center-fed  cylindrical  antenna 
are  given  in  Table  3.  As  can  be  seen  from  this  table,  all  poles  except  the 
first  one  fail  to  meet  criterion  (19) .  The  test  (17)  is  met  by  all  poles 
(thus  (15)  is  violated) ,  so  the  realization  in  the  unmodified  form  (14)  is 
in  no  case  possible. 

Although  the  results  given  in  Table  3  are  somewhat  discouraging,  it  is 
evident  that  the  PR-ness  conditions  are  in  most  cases  only  slightly  vio¬ 
lated.  This  suggests  that,  after  negligible  padding  -  that  is  removing  a 
small  negative  real  part  from  the  function  and  neglecting  it  -  an  approxi¬ 
mate  realization  could  be  possible.  The  results  of  such  an  approach  are 


TABLE  3 


Test  values  for  the  first  layer  poles  given  in  Table  1. 


Test  value 

Pole 

(17) 

(18) 

(19) 

-3 

-3 

-5 

1. 

.1921  x  10 

.9729  x  10 

.1485  x  10 

-3 

-1 

-2 

3. 

.7531  x  10 

.1111  x  10 

-.2945  x  10 

-2 

-1 

-1 

5. 

.1435  x  10 

.3384  x  10 

-.2088  x  10 

-2 

-1 

-1 

7. 

.2187  x  10 

.7004  x  10 

.7048  x  10 

-2 

9. 

.2976  x  10 

.1204 

-.1689 

discussed  in  a  later  part  of  this  report. 

It  is  evident  from  Section  II  that  the  nature  of  the  frequency  character 

istic  of  the  real  part  of  an  immittance  function  Re  f ( jio)  is  extremely 

important  in  the  synthesis  procedure. 

Following  Matthaei  [23],  we  can,  in  general,  classify  immittances  of 

the  form  (6)  in  three  classes:  Class  I  having  Re  f(ju))  .  at  a  =  0;  Class 

min 

II  having  Re  f(jw)  .  at  w  =  and  Class  III  having  Re  f  ( ju>)  .  at  finite 
min  3  J  mm 

values  of  ±  u  (see  Figure  3  (a) ,  (b) ,  and  (c) ,  respectively) . 
o 

As  was  pointed  out  in  Section  II,  an  important  process  in  the  syn¬ 
thesis  of  networks  is  the  procedure  of  removing  a  real  constant1  from  the 
immittance  function,  so  that  the  resulting  remainder  function  f '  (s)  has 
the  property  that 

Re{f ' (jw) }  .  =0.  (20) 


We  allow  the  removal  constant  to  be  either  positive  or  negative,  thus 
implying  the- removal  of  a  resistance  or  conductance  in  the  positive 
case,  or  implying  the  neglecting  of  a  small  negative  quantity  -  i.e. 
padding. 

1 


Lcs  of  three  classes  of  PR  biquadratic 
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Thus,  we  seek  the  frequency  and  the  removal  constant  a  for  which 


mim2 


nin2 


2  2 
m2  -  n2 


s  =  jw 


=  0 


(21) 


min 


Prom  (21)  it  can  be  seen  that 


mim2 


nin2 


am. 


+  an. 


s  =  jto 


=  0 


(22) 


must  have  second-order  roots  at  frequencies  where  (20)  is  satisfied. 
Expressing  (22)  in  terms  of  the  coefficients  of  (6) ,  we  obtain 


D(A  -  aD)u)4  -  (W  -  au)w2  +  F(C  -  aF)  =0  (23) 

where 

W  =  AF  +  CD  -  BE  U  =  2DF  -  E2. 

2 

The  requirement  that  (23)  has  double  order  root  in  u  [47]  leads  to  the 
following  equation  for  a 

(U  -  4D  F  )a  +  (4ADF  +  4CD  F  -  2UW)a  +  (W  -  4ACDF )  =  0  (24) 


The  smaller  of  the  two  roots  of 

X  = 


(24) ,  for  which 


W  -  aU 


2D (A  -  aD)  — 


>  0 


(25) 


is  the  required  solution. 

Having  a,  the  minimum-real-part  frequency  u>  can  be  computed  from 

o 

u  =  /x  .  (26) 

o 


Before  subtracting  a  from  F(s),  we  must  compare  it  with  a^  =  f(0)  and 
aoo  =  f  (°°)  *  The  smallest  value  among  a,  a  ,  and  a^  is  the  sought  minimum 
value  of  Be  f ( jw) . 

The  constant  removal  affects  only  the  numerator  coefficients  of  (6) , 


the  new  values  of  which  are  given  by 


A  =  A'  -  GD 

B  =  B'  -  GE  (27) 

C  =  C'  -  GF 

where  the  old  coefficients  are  primed. 

If  the  a  computed  from  (24)  is  negative,  it  gives  the  value  of  a  shunt 
conductance  (or  series  resistance  -  for  impedance  function)  which  must  be 
added  to  the  non-PR  immittance  f (s)  in  order  to  render  it  a  minimum  PR  func¬ 
tion.  This  procedure  is  called  resistive  padding  [15] . 

As  a  typical  example,  the  real  part  characteristic  of  the  third  unmod¬ 
ified  poler  pair  module  for  the  dipole  antenna  is  given  in  Figure  4.  As  can 
be  seen  from  Figure  4(a)  the  real  part  is  negative  for  low  frequencies. 

This  is  shown  in  detail  in  part  (b)  of  this  figure.  The  minimum  real  part 
frequencies  and  the  resistive  padding  necessary  to  assure  PR-ness  of  the 
unmodified  and  modified  pole-pair  modules  for  the  cylindrical  antenna  are 
given  in  Table  4. 

If  any  of  the  new  coefficients  A,  B,  or  C  in  (27)  happens  to  be  zero, 
the  resulting  function  can  be  easily  synthesized.  Thus,  if  a  =  AVD , 

A  =  0,  and  f(s)  has  the  form  (10).  This  immittance  has  zero  at  infinity 
which  can  be  readily  removed  as  a  single  reactance.  Similarly,  if 
a  =  C'/F,  C  =  0,  and  we  obtain  the  function  of  the  form  (12)  ,  which  has 
zero  at  s  =  0  which  can  be  removed  as  a  pole  of  1/f (s) .  In  both  cases  the 
remaining  function  can  be  easily  realized  with  three  RL  or  RC  elements. 

The  final  networks  for  the  two  cases  with  explicit  expressions  for  the 
element  values,  both  for  f(s)  =  Y(s)  and  f(s)  =  Z(s),  are  given  in  Figures 
5  and  6,  and  Figures  7  and  8,  respectively. 


I 

A 


Figure  5.  Explicit  ladder  networks  for  the  admittance  function  (11)  (zero  at  infinity). 


Figure  6.  Explicit  ladder  networks  for  the  inpedance  function  (11)  (zero  at  infinity) 


Figure  8.  Explicit  ladder  networks  for  the  impedance  function  (13)  (zero  at  zero) . 
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TABLE  4 

Minimum  real  part  frequencies  and  necessary  padding  for 
pole-pair  modules  of  dipole  antenna. 


Pole-pair 

Minimum  real  part 
frequency^- 

Necessary 

padding 

Unmodified 

Modified 

1. 

.0 

2.201 

3  * 

1.598 

4.959 

51.85 

5. 

3.243 

6.507 

32.21 

7. 

4.967 

7.546 

26.44 

9. 

6.708 

8.379 

24.14 

^Normalized. 


The  third  case  possible,  B  =  0,  is  realizable  only  if 

AF  -  CD  =  0  (28) 

(as  a  consequence  of  (19)).  The  resulting  networks  for  the  admittance 
and  impedance  functions  are  shown  in  Figure  9. 

The  situation  is  more  difficult  for  the  case  when  all  the  coefficients 
of  (6)  are  greater  than  zero  but  still  in  some  cases  a  simple  realization 
is  possible.  Thus  if  f(s)  or  its  reciprocal  is  a  Class  I  or  Class  II  PR 
function,  it  can  be  realized  by  a  continued-fraction  expansion  -  a  process 
of  repeated  long  division  and  inversion  applied  to  the  numerator  and  the 
denominator  polynomials  of  the  given  function  and  its  remainders  in  for¬ 
ward  or  reverse  order  [23,22]. 

It  can  be  shown  [23,20]  that  the  biquadratic  function  (6)  can  be 
expanded  in  continued  fraction  only  in  two  cases  when  either 


C 

F 


A 

>  — 
D 


(29) 
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is  guaranteed  by  (33). 
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If 


A 

D 


(and  (33)  holds)  the  following  decomposition  is  possible 

CD  -  AF 


f  ( s)  = 


2 

Ds  +  Es  +  F 


Bs  + 


Ds  +  Es  +  F 


(36) 


(37) 


where  again  the  first  term  is  realizable  as  shown  in  Figure  9.  The  second 
term  has  the  form  (10)  and  is  realizable  as  shown  in  Figures  5  and  6. 

Another  alternative  for  decomposing  a  Class  III  nonminimum  PR  function 
into  sum  of  a  Class  I  and  Class  II  functions  applies  in  some  cases.  It 
takes  the  form 


2 

As  +  Ts  (B  -  T) s  +  C 

2  2 
Ds  +  Es  +  F  Ds  +  Es  +  F 


(38) 


From  (11)  and  (13)  the  range  of  values  of  the  coefficient  T  can  be  deter¬ 
mined  as 


AF  m  BE  -  CD 

—  <  T  <  - - -  . 

E  —  —  E 


(39) 


It  follows  from  (39)  that  this  decomposition  is  possible  only  if 


BE  >_  AF  +  CD  (40) 

-  a  much  stronger  condition  than  (33) . 

It  can  be  shown  that  if  f (s)  fails  to  satisfy  (33)  or  (40) ,  then 
1/f (s)  also  fails  to  satisfy  them. 

A  limited  class  of  biquadratic  minimum  admittances  can  be  realized 
in  a  form  of  five-element  bridge  networks,  as  discovered  by  Kim  and 
VanValkenburg  [19,46]  (see  Figure  10).  The  conditions  which  must  be  satis¬ 
fied  by  the  coefficients  are 


-  29  - 


Figure  10.  Five-e.’ ament  bridge  networks. 
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AF  =  4  BE  =  4  CD 

and 

4  AF  =  4  BE  =  CD 

for  the  networks  (a)  and  (b) ,  respectively.  This  is,  of  course,  a  very 
restrictive  condition.  Thus  these  networks  are  not  of  particular  interest 
in  this  development. 

As  was  mentioned  earlier,  Streable  and  Pearson  [39]  have  shown  that 
the  first  pole  of  the  dipole  antenna  (see  Table  1)  can  be  perturbed  slightly 
(within  a  5%  margin)  so  that  the  modified  pole-pair  admittance  associated 
with  that  pole  is  physically  realizable.  This  admittance  has  the  form  (12) 
and  circuit  (a)  from  Figure  7  can  be  directly  employed.  Unfortunately, 
the  other  pole-pair  modules  of  the  dipole  (and  all  modules,  except  the 
first,  of  the  loop)  fail  to  satisfy  the  PR-ness  condition  (13).  Thus,  a 
resistive  padding  is  necessary.  The  padded  admittances,  however,  has  the 
general  form  (6)  and  the  simple  networks  from  Figures  5  and  7  cannot  be 
used.  It  can  also  be  shown  that  the  decomposition  techniques  discussed 
above  apply  neither  for  the  dipole  nor  the  loop  antenna.  Thus,  the  general 
synthesis  methods  discussed  in  Section  IV  must  be  used. 

It  is  evident  from  Table  3  that  the  PR-ness  test  is  only  slightly 
violated  by  the  third  and  higher  order  modified  pole-pair  modules  for  the 
dipole  antenna.  (Similar  situation  holds  for  the  loop  antenna.)  The 
question  immediately  arises  if  these  pole-pair  admittances  could  be  real¬ 
ized  approximately  by  still  using  the  circuit  (a)  from  Figure  7  and 
neglecting  small  negative  elements  which  could  arise.  Such  attempt  has 
been  made  and  the  network  shown  in  Figure  11  was  obtained.  The  transient 
response  of  this  network  for  the  Gaussian  pulse  defined  in  Figure  12  was 
evaluated  by  SCEPTRE  [9]  circuit  analysis  program  and  compared  with  the 


Figure  11.  Approximate  equivalent  network  for  dipole  antenna. 


TWTD 


Figure  13.  Transient  response  of  the  lumped  circuit  (Figure  11) 
and  the  result  of  the  TWTD  program  analysis. 


response  obtained  from  the  TWTD  [25]  thin  wire  analysis  program.  The 
results  are  shown  in  Figure  13.  It  is  evident  that,  although  the  main 


features  are  preserved  in  the  circuit  response,  the  agreement  between  the 
two  curves  is  not  satisfactory. 

Another  approximate  equivalent  network  for  the  dipole  antenna  results 
from  the  approach  suggested  by  Baum  and  Singaraju  [6] .  They  approximated 
the  natural  modes  of  the  dipole  by  sinusoidal  functions  and  derived  analy¬ 
tically  the  following  expression1  for  the  residues  (real  and  equal  for  all 
poles) 


o 


where  ft  is  the  thickness  factor  and  Z^  is  the  intrinsic  impedance  of  free 
space.  For  the  dipole  from  Table  1  a  1  .001.  For  this  residue  and  the 
poles  listed  in  Table  1  the  PR-ness  condition  (14)  is  always  satisfied  and 
the  input  admittance  is  realizable  in  the  form  (10) .  Using  the  circuit 
(a)  from  Figure  5  we  obtain  the  network  shown  in  Figure  14.  One  defic¬ 
iency  of  this  network  is  apparent:  its  admittance  is  finite  at  zero  fre¬ 
quency,  which  is  inconsistent  with  the  structure  under  consideration. 

The  results  of  the  SCEPTRE  and  the  TWTD  analyses  are  shown  in  Figure 
15.  The  same  Gaussian  pulse  as  before  was  used. 

Again,  the  agreement  between  these  curves  is  not  very  good,  but  the 
last  approach,  as  more  systematic  than  the  previous  ad  hoc  procedure,  is, 
of  course,  more  satisfactory. 


Normalized  to  rc/L. 


i 


Equivalent  network  for  the  dipole  antenna  resulting 
from  the  approximation  (41) . 


Section  IV 


SYNTHESIS  OF  MINIMUM  BIQUADRATIC  FUNCTIONS  -  GENERAL 

METHODS 


It  is  shown  in  Section  III  that  the  equivalent  network  synthesis  on 
pole-pair  at  a  time  basis  leads  in  general  to  nonrealizable  biquadratic 
modules  for  the  dipole  and  loop  wire  antennas.  The  real  part  characteristic 
shown  in  Figure  4  is  typical  for  these  modules.  It  is  also  shown  that  the 
PR  character  of  pole-pair  admittances  can  be  restored  by  negligible  resis¬ 
tive  padding.  It  is  conceivable  that  this  also  holds  for  other  (at  least 
high-Q)  structures. 

After  resistive  padding  we  obtain  a  minimum  biquadratic  function.  Thus, 
the  function  (6)  for  which  the  condition 

(/ftF  -  /CD)  =  BE  (42) 

holds.  In  general,  this  will  be  a  Class  III  function  (see  Figure  3) ,  the 
real  part  of  which  is  zero  at  some  finite  frequency  (if  it  happens  to 
be  a  Class  I  or  Class  II  function  it  can  be  easily  synthesized  by  the 
methods  of  Section  III).  From  (8)  and  (42)  it  can  be  shown  that 
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4.1  BRUNE  METHOD 

Intrinsic  to  this  method  is  the  use  of  a  perfect  transformer.  Thus,  it 
is  rather  impractical  and  is  seldom  used,  except  as  a  last  resort.  On  the 
other  hand,  the  Brune  network  contains  the  minimum  number  of  elements  possi¬ 
ble  and  thus  can  be  convenient  for  analysis  and  modeling  purposes. 

At  the  frequency  (i>o  the  whole  impedance  Z(jw)  =  1/Y(jw)  is  purely 
reactive 


z(jw)  =  jx«oo>  =  yBTuTT 


and  can  be  represented  by  a  single  inductance 


(47) 


(48) 


S' 
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For  the  case  A,  using  (43)  and  (46a)  we  obtain 


/DE 

L1  /bc 


(49) 


Although  is  negative  we  proceed  to  remove  the  impedance  sL^  from  Z(s): 


^  ^  >  .  Ds  +  Es  +  F  /DE 

Zi  (S)  =  Z  (S)  -  SL1  =  2 -  +  /ic  3  ' 

As  +  Bs  +  C 


(50) 


-  2  2 
We  know  that  Z(s)  has  a  zero  at  s  =  ± jui  thus  the  term  (s  +  w  )  can  be 

o  o 

factored  out  in  (50)  to  give 


~  ,  ,  ,2  2S  /  BC 

Z  (s)  =  (S  +  W  )  - r— 

1  0.2 


.  /DE  ±  n  /5f 

A  /—  s  +  A  /  — 
/  AC 


As  +  BS  +  C 


(51) 


The  zero  of  Z^  can  be  removed  as  a  pole  from  Y^(s)  =  l/Z^Cs) 


2a  s 

?1(S|  ’  T  P  2  +  Vs) 

s  +  0> 


(52) 


where  the  first  term  represents  the  admittance  of  a  series  I^C^  circuit 
(Figure  16) .  The  residue  aQ  can  be  easily  computed  as 


.  B  2  B  /CF 

2a  =  —  o)  =^/rr- 
o  F  o  F  v  AD 


(53) 


The  values  of  L.  and  C  are  given  by 
2  o 


1  _  F  1  _  F  /AD 

L2  ~  2a  “  B  2  B  /  CF 
o  w 

o 


(54) 


C  = 
o 


2. 

“o  L2 


B 

F 


(55) 


From  (52)  and  (53)  the  impedance  Z^  =  1/Y ^  is  obtained  as 
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z2  (S) 


which  can  be  realized  as  a  series  L^R  circuit,  where 


(56) 


L 


3 


(57) 


and 


R  = 


F 

C  * 


(58) 


The  resulting  network  is  shown  in  Figure  16.  The  negative  inductance  can 
be  realized  with  the  help  of  a  perfect  transformer.  The  network  shown  in 
Figure  17  results. 

By  development  analogous  to  the  above,  it  can  be  shown  that  this  network 
is  also  valid  for  the  case  E. 

As  an  example,  the  Brune  equivalent  network  for  the  dipole  antenna  is 
shown  in  Figure  18.  (The  amount  of  resistive  padding  used  is  given  in  Table 
4.  For  the  first  pole-pair  admittance,  which  is  a  Class  I  function,  the  net¬ 
work  (a)  from  Figure  7  was  used.)  The  transient  response  of  this  network 
for  the  Gaussian  pulse  (Figure  12)  is  compared  with  the  result  of  the  TWTD 
analysis  in  Figure  19.  Except  for  the  early  time,  an  excellent  agreement  is 
observed . 


4.2  BOTT-DUFFIN  METHOD 

This  is  the  most  general  transformerless  synthesis  method.  It's  main 
disadvantage  in  general  application  is  the  exponential  growth  of  the  number 
of  elements  required  with  the  number  of  cycles  needed  to  reduce  the  given 
function  to  a  constant  remainder.  This  difficulty  is  insignificant  in  the 
present  application  of  synthesis  of  biquadratic  functions.  Therefore  it  is 


Figure  18.  Brune  equivalent  network  for  dipole  antenna 


Transient  response  of  the  lumped  circuit  (Figure  18)  and 
the  result  of  TVJTD  program  analysis. 


considered  in  detail  here. 


In  the  Bott-Duffin  procedure,  one  makes  use  of  the  Richards  function 


kY(s)  -  sY(k) 
kY  (k)  -  sY(s) 


(59) 


in  which  k  is  a  real  positive  constant.  It  may  be  shown  that  R(s)  is  PR  if 
Y (s)  is  PR,  and  that  the  order  of  it  (defined  as  the  total  number  of  its 
poles  or  zeros)  is  not  higher  than  that  of  Y(s). 

The  equation  (59)  can  be  solved  for  Y(s)  to  give 


Y(s) 


Vs)  +  ir  +  22(s)  +  SLo 

o 


(60) 


where 


Z  (s) 


Z2(s) 


ft(s) 
Y  (k) 


1 

Y (k)  R(s) 


(61) 


(62) 


r  - 
o  k 


(63) 


Lo  kY(k) 


(64) 


This  corresponds  to  the  network  shown  in  Figure  20.  Note  that  the  PR-ness 

of  Y (s)  assures  that  the  values  C  and  L  are  positive  and  real. 

o  o 

It  is  easy  to  see  that  the  network  from  Figure  20  is  a  balanced  bridge 
with  the  points  A  and  B  on  the  same  potential.  Thus,  any  impedance  (even  a 
short  circuit)  can  be  placed  in  the  detector  arm  A-B  without  affecting  the 
input  admittance  of  the  circuit.  This  can  be  exploited  to  reduce  by  one 
the  number  of  reactive  elements  necessary  [27,32,37], 


d 
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At  the  frequency  s  =  from  (60)  ,  we  have 


Y(jw  )  =  jB(ui  )  =  - - - 

°  °  2(jw)+-i~ 

1  o  noi  C 
o  o 


Z2(^0)  +  ^0Lq 


The  special  cases  A  and  B  referred  to  in  the  preceding  section  are 
treated  separately. 

Case  A 

Because  Bfu^)  >  0  for  this  case,  the  only  possibility  in  satisfying 

(65)  occurs  when  Z  ( ju)  )  =  0  and  2_(jio  )  =  From  (61)  and  (62)  we  see  that 

1  o  2  o 

this  requires  that  R(s)  have  a  zero  at  s  =  ju  .  This  in  turn  requires  that 

o 

the  numerator  of  (59)  be  equal  to  zero  at  this  frequency: 


kY  ( j )  -  ji^Yfk)  =  0 


w  Y(k)  -  kB(u>  )  =  0. 
o  o 


From  (67)  with  the  help  of  (46a)  the  value  of  Cq  can  be  computed  as 

_  Y  (k)  _  S<V  _  Ab 

o  ~  k  0)  /  EF  ' 


In  order  to  determine  the  value  of  L  ,  eq.  (67)  must  be  solved  for  the  real 

o 

positive  root  k: 


Ak  +  Bk  +  C  ,  /AB 

-1 - k  /—  -  0. 

Dk  +  Ek  +  F 


With  the  help  of  (44a)  this  equation  leads  to 


,  3  /CE  ,  2  /CF  ,  C  /EF 
k  "  /  IFk  +  /  AD  k  ”  D  /  AB 


0. 


(70) 
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Although  (70)  is  a  third  order  equation,  the  difficulty  is  an  apparent  one 

only,  because  we  know  the  two  of  the  roots:  k  =  ±jwo.  Thus,  dividing  (70) 
2  2 

by  (k  +  to  )  we  find 


k  = 


(71) 


From  this  equation  and  from  (68)  it  follows  that 


Y(k) 


(72) 


and 


_  1  _  D  /5F 

o  kY  (k)  C  /  AE 


(73) 


The  impedance  Z^s)  has  a  zero  at  s  =  jo)Q  and  Z2(s)  has  a  pole  at  this 
frequency.  The  poles  of  1/Z^(s)  and  Z^(s)  can  be  removed  by  the  following 
expansions 


1 

Z1  (s) 


2a^s 


2 

s  +  0) 


2 


o 


(74) 


2a  s 

Z2(S)  =  1 - 2  +  V 

s  +  w 


(75) 


The  residues  a^  and  a2  determine  the  values  of  the  series-connected  and 
and  the  parallel-connected  C2  and  L2  (see  Figure  21) . 

In  order  to  compute  the  residues  a^  and  a2  we  note  that  the  numerator 
and  denominator  of  the  Richards  function  must  have  the  common  factor  (s-k) 
[49].  After  dividing  it  out,  with  the  help  of  (44a),  (71),  and  (72), 


we  get 


R(s) 


(76) 
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m  =  b-%/df  +  eVac 


Figure  21.  Bott-Duffin  network  for  the  Case  A. 
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and 


where 


1 

^(s) 


_M _ 

D^DF 
2  ^ 

S  +  ID 

O 


M  =  B/dF  +  e/aC  . 

Comparing  (77)  and  (78)  with  (74)  and  (75)  we  obtain 


(77) 


(78) 


(79) 


C 


1 


(80) 


1  d/df 

1=2ai“  M 


c 


2 


L 


2 


(81) 


(82) 


(83) 


R 


1 


(84) 


R 


2 


(85) 


The  final  network  for  the  case  A  is  shown  in  Figure  21. 

Case  B 

In  this  case  8(<dq)  <  0  and  (44b)  and  (46b)  apply.  From  (65)  we  see 
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r 


that  1/Z^  ( jto  )  =  0  and  Z^juM  =  0  which,  by  virtue  of  (61),  (62),  and  (59), 


leads  to  the  condition 


kY  (k)  -  Y  ( jw  )  =  0 

o  o 


kY  (k)  +  w  8(fa>  )  =  0 
o  o 


From  (65)  and  (87) ,  using  (43)  and  (46b)  we  obtain 


T  =  =  __=1 _  =  /“ 

o  kY  (k)  fa)  8  (w  )  /  BC 

o  o 


In  order  to  determine  the  value  of  C  equation  (87)  must  be  solved  for  k. 

o 

With  the  substitutions  (43)  and  (46b)  we  obtain 


+  Bk  +  C  /BC 
_  ~  /  DE 


Dk  +  Ek  +  F 


which,  with  the  help  of  (44b)  leads  to  the  equation 


k3  .  rn.  k2  +  M.  *  .  £  /2£  .  0. 

/  AE  /AD  A  /  DE 


2  2 

Dividing  (90)  by  (k  +  u)q  )  ,  as  in  case  A,  we  get 


k  =  /!!• . 

/  AE 


Using  (88)  we  also  have 


Y<k>  =  M 


A  /CE 


F  /  BD 


Y  (k) 
k 


(93) 
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The  Richards  function  can  be  expressed  as 


R(s)  = 


^£s2  +(*  + 

AF _ \A  /  DF  AJ 


(94) 


2 

s  + 


'CF 

AD 


The  poles  of  Z^(s)  and  l/Z^s)  can  be  removed  similarly  as  in  case  A 


2a  s 

VS)  =~2 - 2 

s  +  U) 


+  R, 


(95) 


2a2s 


22 (s)  s  +  u 


2  R„ 


(96) 


Using  (61) ,  (62) ,  and  (94)  we  obtain 


Vs)  -  2 


M 

A^AC  S 


D 

+  T  * 


s  +  0) 


2  A 


(97) 


and 


M _  /C 

1  _  df/  a  s  c 

+  T-l  ' 


Z,(s)  2  2  ■  F 

2  s  +  w 

o 


(98) 


where  M  is  defined  as  in  case  A.  Comparison  with  (95)  and  (96)  gives 


1  _  a/Kc 

C1  2a:  =  M 


(99) 


-  2h.  _  M_ 

1  2  AC/  F  ' 

0) 

o 


(100) 


2a„ 


C2  = 


2  f/dF 


(101) 


Figure  22.  Bott-Duffin  network  for  the  Case  B. 
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L 


2 


1  _  DF  /K 

2a2  M  /  C  ' 


(102) 


*1  -  |  '  (103) 

and 

»2*§-  <1MI 

The  final  network  for  the  case  B  is  shown  in  Figure  22. 

Modified  networks 

As  was  mentioned  previously,  in  both  cases  A  and  B  one  element  can  be 

spared  by  placing  appropriate  reactance  in  the  detect  detector  arm  A-B  and 

performing  delta-T  or  T-delta  transformation.  It  has  been  shown  by 

Balabanian  and  Cahn  [2]  that  in  order  for  the  T-delta  transformation  to  be 

physically  realizable  the  impedance  Zb(see  Figure  23)  must  be  inverse  to 

Z  and  Z  must  be  inverse  to  the  sum  of  the  first  two.  similar  conditions 
c  a 

must  hold  for  the  delta-T  transformation.  It  is  easy  to  see  that  the  net¬ 
works  shown  in  Figure  24  satisfy  the  above  conditions  if  the  resonance  fre¬ 
quencies  of  the  circuits  formed  by  and  are  equal.  Thus  if 


L  C 
a  a 


L  C. 
c  b 


(105) 


The  relations  between  the  elements  in  the  T  and  Delta  networks  are  also 
given  in  Figure  24  [11]. 

Consider  the  case  A  first.  From  the  above  considerations  it  follows 
that  a  capacitor  should  be  placed  between  the  points  A  and  B  in  Figure 
21,  the  value  of  which  is  given  by  condition  (105)  as 


C 


b 


A  /CE 

f/  bd 


(106) 


Figure  23.  Illustration  of  the  T-delta  transformation. 


Lq 


Figure  24.  T-delta  transformation  used  to  simplify  the  Bott-Duffin  networks 
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After  replacing  the  T  which  appears  to  the  right  (Figure  21)  by  delta,  the 
network  shown  in  Figure  25(a)  results.  The  equivalent  network  given  in 
Figure  25(b)  is  obtained  when  the  delta  which  appears  at  the  bottom  in  Fig¬ 
ure  21  is  replaced  by  T. 

For  the  case  B,  by  analogous  procedure,  by  placing  an  inductor 


L 

c 


D  /BF 
C/  AE 


(107) 


in  the  detector  arm  of  the  network  in  Figure  22,  we  obtain  the  two  equiva¬ 
lent  modified  Bott-Duffin  networks  given  in  Figure  26. 

It  should  be  pointed  out  that,  although  the  two  networks  given  in 
Figure  25  are  fully  equivalent,  one  of  them  may  have  more  desirable  element 
values  from  the  realization  point  of  view.  The  same  holds  for  the  networks 
in  Figure  26. 

As  cam  be  seen  in  Figures  25  and  26,  the  Bott-Duffin  networks  even  in 
the  reduced  form,  are  more  complicated  than  the  Brune  network  (Figure  17) . 
Apparently  the  price  we  pay  to  avoid  transformers  is  extra  elements. 

Simplified  Bott-Duffin  Networks 

It  is  of  interest  to  note  that  all  pole-pair  modules  for  the  dipole 
and  loop  wire  antennas  qualify  for  the  case  A,  i.e.  at  the  zero  real  part 
frequency  they  can  be  approximately  represented  by  a  single  capacitance 
Cq  given  by  (68) . 

Real  part  characteristic  of  a  typical  modified  pole-pair  module, 
before  and  after  padding,  is  shown  in  Figure  27.  From  Figures  27  and  21 
it  is  evident  that  R^  =  1/G.  If  the  padding  G  is  really  negligible,  R2 
is  a  large  resistance.  Consequently,  the  whole  branch  in  series  with  R2 


contributes  insignificantly  to  the  total  admittance  and  can  be  neglected. 


R  =  R/BD  +  Ev€e 


M  -  BVf)F+Ev£C 

(a) 


Figure  25.  Modified  Bott-Duffin  networks  for  the  Case  A. 


c2=24Cv'S§ 

^  m2  AF 

>  CDV  BF 
Q  =  AVcE  +  BTBD 


M  =  Bv'DF  +  EVKC 


Figure  25  (Continued) . 


r 


Figure  26.  Modified  Bott-Duffin  networks  for  the  Case  B. 


M  »  bTdF  +  Ev^C 

(b) 


Figure  26  (Continued) . 
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This  simplification  was  first  discovered  by  Streable  and  Pearson  [39]. 

It  is  conceivable  that  for  some  structures,  made  of  a  non-perfectly  con¬ 
ducting  material,  the  real  part  characteristic  of  pole-pair  admittances  has 
the  form  shown  in  Figure  28.  Then,  if  the  case  B  applies,  the  resistor 
R^  =  1/G  in  Figure  22  is  large  and  the  whole  branch  in  series  with  R^  can  be 
neglected. 

The  simplified  Bott-Duffin  networks  for  both  cases  A  and  B  are  shown 
in  Figure  29. 

The  admittance  of  the  simplified  Bott-Duffin  network  for  case  A  (see 
Figure  29  (a) )  is  given  by 


This  is  a  third  order  function  which  approximates  the  original  second  order 
admittance  (6).  Note  that  besides  the  original  poles,  an  additional  real 
pole  (which  happens  to  be  close  to  the  origin  for  thin  wire  dipole  and  loop 
antennas)  is  introduced.  Also  a  zero  at  the  origin  is  added.  It  is  evident 
that  for  higher  frequencies  the  added  zero  and  pole  will  approximately  can¬ 
cel  and  this  function  will  differ  insignificantly  from  the  original  one. 

The  equivalent  network  for  the  dipole  antenna  using  the  simplified  Bott- 
Duffin  module  from  Figure  29  (a)  is  shown  in  Figure  30. 1 

The  transient  response  of  this  network  for  the  Gaussian  pulse  (Figure  11) 
is  compared  in  Figure  31  with  the  result  of  TWTD  analysis.  The  agreement  is 
excellent  except  for  the  early  time.  Corresponding  results  for  the  loop 
antenna  are  shown  in  Figures  31A  and  31B. 

■'"The  first  pole-pair  module  was  realized  by  the  method  given  in  Section  III. 


Simplified 


Transient  response  of  the  lumped  circuit  (Figure  30) 
and  the  result  of  TWTD  program  analysis. 


8.94 


Figure  31B.  Transient  response  of  the  lumped  circuit  (Figure  31A) 
and  the  result  of  the  TWTD  program  analysis. 
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4.3  DARLINGTON  METHOD 

In  this  method  we  seek  the  realization  of  a  given  function 


Z(s) 


mi  +  ni 
m2  +  n2 


(108) 


in  a  form  of  an  input  impedance  of  a  lossless  two-port  terminated  in  a  single 
resistor  (Figure  32).  Thus,  Z(s)  can  be  written  in  the  form 


Z(s)  = 


Rz^1  +  Izl 


1  + 


Izl 


Rz 


R  +  Z 


=  z 


11 


22 


11  z 
1  + 


22 

R 


(109) 


where  the  z.  ,'s  are  the  open-circuit  admittance  functions. 

13 

By  comparison  with  (109),  (108)  can  be  rewritten  in  two  ways  as  follows: 


11+  n  /m 

Z(s)  =  —  -----  ;  (case  A)  (110) 

n2  1  +  m2/n2 


and 


n  1  +  m  /n 

Z(s)  =  —  — — (case  B)  .  (Ill) 

m2  1  + 


Thus,  we  can  make  the  following  identifications. 


Case  A 
m. 


Z11  n. 


Case  B 
n. 


'll  m„ 


(a) 


'22  _  m2 
R  n^ 


'22  n2 
R  m_ 


(b)  (112) 


'12 

R 


/l(-s2) 


'12 

R 


/n(-S2) 


(c) 
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with 

N(-s2)  =  m^ir^  -  n^n2  .  (113) 

It  is  easy  to  see  that  z ^  and  z ^ ,  being  the  ratios  of  odd  and  even 

parts  of  Hurwitz  polynomials,  are  realizable  reactance  functions.  In  order 

to  satisfy  the  conditions  of  realizability  of  a  lossless  two-port  we  must 

2  2 

require  that  z ^  be  a  rational  function  (N(-s  )  or  -N(-s  )  must  be  a  complete 
square)  and  the  residue  condition 


k 


llk22 


(114) 


be  satisfied  at  all  the  poles  [1,16]. 

For  all  finite  poles  s^  =  in  the  case  A,  this  condition  takes  the 

form 


(115) 


Thus,  the  residue  condition  is  satisfied  with  the  equality  sign  at  all  finite 
poles  (the  poles  are  compact) .  Similar  conclusion  holds  for  case  B. 

The  set  of  three  rational  functions  z  ,  z22'  an<^  Z12  ^or  t*ie 

residue  condition  (114)  is  satisfied  at  each  pole  can  be  synthesized  by  the 
Cauer  realization  procedure  [1,49].  In  this  method  the  set  of  z  parameters 


(112)  is  expanded  in  partial  fractions  as  follows 

r  l  r 


'11 


22 


12 


v  1«»  • 

-ii  I  * 

s 

I 

(0)  | 


Rk 


22 


/R  k 


(0)| 

12.  I  + 


(®) 

ku  . 


(“) 

**22  ’s 


u 


..  (1) 
2kll  3 


l+l 


/r  k 


12 


s .  + 


I  I 


2 

2 

s 

+ 

2Rk  (Da  1 

22 

2 

2  1 

s 

H 

3 

+ 

/r 

k 

12  1 

2 

2 

s 

+  "l  J 

(116) 
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where  only  one  pair  of  finite  j-axis  poles  has  been  explicitly  shown.  Next, 
the  parameters  in  each  of  the  brackets  are  realized  as  a  two-port  and  the  com¬ 
plete  realization  is  obtained  as  a  series  connection  of  these  components. 

If  the  poles  are  compact,  the  component  two-port  can  be  synthesized  by 
the  network  shown  in  Figure  33.  The  open-circuit  impedance  functions  for 
this  network  are 

2 

z^  =  Z,  z22  =  n  Z,  and  z12  =  nZ  .  (117) 

The  parameters  in  any  one  of  the  brackets  in  eg.  (116)  can  be  represented  as 


and 


211  =  klig(E)' 


Z22  =  ^22g{5)' 


Z12  =  ^  ki2g(s)' 


2  2 

where  g(s)  represents  1/s,  s,  or  2s/(s  +  ) . 


Comparison  of  (117)  and  (118)  gives 


(118) 


Z  =  kj^g(s)  ,  (119) 

and  , 

—  k?-> 

n  =  V*R  T~~  •  (120) 

k12 


It  is  evident  from  (120)  that  the  ideal  transformer  associated  with  the 
given  pole  can  be  avoided  by  choosing 


R  = 


(121) 


Returning  to  the  biquadratic  function  (6)  we  have 

2 

=  Ds  +  F,  n^  =  Es, 

2 

m2  =  As  +  C,  n2  =  Bs,  (122) 

2  4  2 

N ( -S  )  =  ADs  +  (AF  +  CD  -  BE)s  +  CF. 


and 


Making  use  of  the  minimum  condition  (42) ,  the  last  expression  can  be  written  in 
the  form 

N (— s  )  =  (►'AD  s  +  /CF  )  .  (123) 

2 

Thus,  N(-s  )  is  a  complete  square  and  the  case  A  applies. 


Substituting  (122)  into  (112)  we  obtain 


FI  D 

Z11  B  s  +  B  S 


_  C  1  ,  _  A 
*22  ’  R  B  I  +  B  B  S 


(b)  (124) 


for  the  pole  at  zero  and 


for  the  pole  at  infinity. 

The  Cauer  realization  is  shown  in  Figure  34.  The  presence  of  perfect 
transformers  makes  this  network  highly  undesirable  for  practical  realization. 
The  upper  ideal  transformer  can  be  however  replaced  by  a  perfect  transformer 
according  to  the  equivalence  shown  in  Figure  35.  The  lower  ideal  transformer 
can  be  completely  eliminated  if  we  choose  R  according  to  (121).  Thus, 

R  =  F/C.  The  final  network  is  shown  in  Figure  36.  Being  identical  with  the 


Brune  network  derived  previously,  it  suffers,  of  course,  from  the  same 
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deficiencies. 

Hazony  [16]  described  a  modification  of  the  above  procedure  leading  to 
no  transformers  at  the  expense  of  more  RLC  elements.  This  approach  is  dis¬ 
cussed  below. 

Consider  the  even  part  of  a  minimum  biquadratic  impedance  function 


Z(s)  =  1/Y (s) 


Ev{Z(s) } 


2  2  2 
AD  (s  +  0)  ) 

_ _ o _ 

2  2  2  2 
(As  +  C)  -  B  s 


(125) 


where  either  (44a)  or  (44b)  holds  and  ca  is  given  by  (43)  . 

o 

2  2  2  2 

By  multiplying  (125)  by  the  surplus  factor  (k  -  s  )/(k  -  s  )  it  can 

be  decomposed  as  shown  below 


2  2  2 
k  AD  (s  +  w  ) 
o 


Ev{Z  (s)  }  =  - j - 2 - 1~ 2 - ' 

[  (Ak  +  B)s  +  kC]  -  S  [As  +  C  +  Bk] ' 


2,2  22 
-  s  AD  (s  +  u  ) 
o 


[ (Ak  +  B)s2  +  kC]2  -  S2[AS2  +  C  +  Bk]2 


(126) 


or 


Ev{Z (s) }  =  Ev(Z1(s)}  +  Ev{Z2(s)}  =  — y 


m2  "  n2 


2  ' 


(127) 


where  Z^(s)  and  Z2(s)  can  be  synthesized  by  Darlington's  procedure.  The 
positive  constant  k  can  be  chosen  so  that  no  transformers  result. 
Consider  Z^(s)  first.  From  (112)  we  have 


"22 


m 


‘2  _  (Ak  +  B)s  +  kC 
2 

s(As  +  Bs  +  C) 


(128) 


m 

o 


k/AD( 
s  (As2 


2 

s 


+ 


+  C  +  Bk) 


(129) 
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There  are  two  distinct  cases  which  will  lead  to  no  transformers:  (1)  the  zero 
of  z  will  cancel  its  pole  or  (2)  z  has  same  zeros  and  poles  as  z  . 

X.C.  ££ 

For  the  first  case  we  have 


C  +  Bk  2 

— - -  =  to 

A  O 


or 


k  = 


/AF  -  /CD  /C 


D  ' 


which  is  positive  only  when  (44a)  holds  (case  A) .  Thus, 


k  =  m 

/  BD 


(130) 


In  the  second  case  above  we  have 


kC 


Ak  +  B  o' 


(131) 


or 


k  = 


✓CD  -  /AF  v  A 


which  is  positive  only  when  (441  holds  (case  B) .  Thus, 


k-/I- 


(132) 


In  both  cases  the  remaining  impedance  z  can  be  found  from  the  residue 
condition  (the  residues  must  be  compact) . 

Consider  the  case  A  first.  Substituting  (130)  for  k  in  (128)  and  (129) 
we  obtain  (after  partial  fraction  expansion) 


C  /EF  M 

Z22  =  F/  AB  S  aTdf 

K  ~  S  +  2  2 

i  s  +  to 


(133) 
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/CE 

Z12  /aB 


and,  from  the  residue  condition, 


with 

M  =  Bv/DF  +  e/aC  . 


(134) 


(135) 


This  set  of  parameters  can  be  now  easily  realized  by  the  Cauer  method 
described  previously.  For  the  pole  at  zero  we  eliminate  the  ideal  transformer 
by  choosing  the  value  of  R 2  as  prescribed  by  (121).  Thus,  R2  =  F/C.  The 
pole  at  ±j“Q  is  a  private  pole  [47]  of  z22  and  can  be  simply  realized  as  a 
parallel  LC  network  connected  in  series  with  the  output.  The  resulting  net¬ 
work  is  shown  in  Figure  37. 

Consider  the  case  B  next.  From  (128) ,  (129) ,  and  (132)  we  obtain 


c/bf  mb  /c 

R  S  RA/  E 

s  2  E 

s  +  av'Se 


and 


(136) 


s 


2 

s 


MB  /¥_ 

RA/  E 
- 

+  a/ae 


(137) 


with 


R  =  CVAE  +  B/iF. 
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= 


_  R 


FJBF 


B  _  F 

R2“T 


C2 


raVce 

MBF 


^-“54 


R  -  cVAE  +  B^/BF  M  *  By/?)?  +  E-/AC 


Figure  37.  Darlington  procedure: 


realization  of  Z^(s) 


Case  A 
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From  the  residue  condition  we  get 


z 


11 


f/If 


MBF 

S  RA/(5e 
2  _R_ 

3  +  A*5£ 


(138) 


The  choice  =  F/C  leads  to  a  transformerless  realization  shown  in  Figure 
38.  This  network  is  equivalent  to  network  shown  in  Figure  39. 

For  the  impedance  Z  (s)  ,  from  (126),  (127),  and  (112),  we  have 


Z22  n2  _  s (As2  +  C  +  Bk) 
R1  m2  (Ak  +  B)s2  +  kC 


and 


n  s/AD  ( s  +  co  ) 
o  o 


'22 


12 


and 


'11 


m2  (Ak 

2 

+  B)  s  +  kC 

For  the 

case  A  we  get 

a/bd  , 

MB 

s  ;t/ac 

“  Q 

2  ,  C/CE  ' 

3  Q 

d/ab 

MB 

S  2/CD 

Q 

“  Q 

2  C/CE  ' 

3  +  Q 

DV'BD 

MBD  /c 

S  2  /  A 

Q 

“  Q 

2  ,  C/CE 

3  Q 

(139) 


(140) 


(141) 


(142) 


(143) 


where 


Q  =  A/CE  +  B/BD  . 


The  choice  R^ 


=  D/A  leads  to  a  transformerless  realization  shown  in  Figure  40. 
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This  network  is  equivalent  to  the  network  shown  in  Figure  41. 


For  the  case  B  we  get 


A  /BE  .  D/AC 
S  b/cD  ~ - 2  ' 

S  +0) 

o 


2  =  S 

11  V 


With  =  D/A  this  set  of  parameters  has  the  realization  shown  in  Figure  42. 

If  we  now  connect  in  series  the  networks  for  Z^  and  Z^  for  the  cases  A 
and  B,  the  networks  in  Figure  43  result. 

It  can  be  shown  by  using  the  equivalent  networks  from  Figures  38  and  41, 
that  the  networks  obtained  are  the  Bott-Duffin  networks  (Figures  21,  22) 
except  that  the  detector  arms  are  short-circuited. 


4.4  MIYATA  PROCEDURE 

In  this  method  the  even  part  of  the  given  function  Z(s)  is  decomposed 


as  follows 


Ev{Z (s) }  = 


N(-S2) 


„  ,  2.P 

T  a  (-s  ) 


2  2  2  2 
m2  n2  m2  n2 


,  2s  /  2\n 

a  a,  (-s  )  a  (-s  ) 

.  o  +  _ +  +  _2 _ 

2  2  2  2  2  2 
m2  n2  m2  n2  m2  “  n2 


Ev(Z  }  +  Ev{Z, }  +  ...  +  Ev{Z  } 
o  1  n 


Figure  43.  Final  form  of  Darlington  networks  for  Cases  A£a),  and  B(b) 
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r  M  S 

0  FVBF 


R  *  C>/AE  +  BVBF 


r  RAv^CE 
MBF 

Vaf  -  v'cd  =  -Tbe 


(b) 


Figure  43  (Continued) . 
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If  the  coefficients  a^  are  positive,  each  of  the  even  parts  in  (147)  will 
correspond  to  an  impedance  which  can  be  realized  as  a  network  containing  n 
reactive  elements  plus  a  resistance. 

Even  if  some  of  the  coefficients  a^  in  the  original  function  are  negative 
they  can  be  made  positive  by  multiplying  the  numerator  and  the  denominator 

of  Z (s)  by  suitable  surplus  factors  [l].  This  procedure  fails  however  if 

2  2  2  2 
N(-s  )  contains  factors  of  the  form  (s  +  u  )  .  A  look  at  eg.  (125)  reveals 

o 

that  our  function  falls  into  this  very  class.  Thus,  the  Miyata  procedure 
is  not  applicable  to  minimum  functions. 
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Section  V 


PRACTICAL  MATTERS 

5.1  CAPACITIVE  ADJUSTMENT 

It  is  of  interest  to  point  out  that  all  the  SEM  derived  equivalent  net¬ 
works  given  in  Section  IV  fail  to  reproduce  the  early  time  response  of  the 
structure,  as  computed  by  TWTD  program.  It  was  suggested  by  Streable  et.  al. 
[40]  that  this  departure  of  results  is  attributable  to  the  quasistatic  influ¬ 
ence  of  the  reactance  which  one  neglects  when  the  summation  (1)  is  truncated 
to  a  finite  number  of  elements  for  realization.  This  asymptotic  reactance  of 
a  given  structure  must  be  computed  apart  from  the  SEM  representation  and  the 
reactance  of  the  truncated  equivalent  circuit  must  be  adjusted  with  additional 
element (s)  so  as  to  match  that  of  the  object  in  the  asymptotic  limit. 

As  an  example  the  static  capacitance  of  the  simplified  Bott-Duffin  net¬ 
work  for  a  dipole  antenna  (Figure  30)  was  adjusted  by  adding  a  70  pF  capaci¬ 
tance  in  shunt.  A  remarkable  improvement  in  the  early  time  response  of 
the  network  was  observed  (Figure  44) .  This  particular  value  of  C  was  found 
by  a  cut  and  try  method  but  systematic  means  of  finding  it  are  conceivable. 

For  example,  C could  be  determined  as  a  difference  between  the  static 
capacitance  of  the  structure  determined  by  some  auxiliary  means  -  for 
example,  a  method  of  moments  solution  for  the  static  capacitance  -  and  the 
static  capacitance  of  the  equivalent  network  synthesized  from  the  SEM  data. 
Unfortunately,  it  is  believed  that  this  method  cannot  be  employed  in  the 
case  of  approximate,  resistively  padded  networks. 

5.2  PARASITIC  EFFECTS 

All  networks  considered  thus  far  comprise  only  ideal  RLC  elements. 
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TWTD 


Figure  44.  Early  time  response  of  the  equivalent  network 
(Figure  30) .  (a)  Without  correction,  and  (b) 


However,  in  the  practical  realization  each  element  introduces  some  parasitic 
elements  which  have  deteriorating  influence  on  the  network  performance.  There 
fore  the  impact  of  the  parasitic  effects  should  be  carefully  considered, 
particularly  if  the  network  has  to  operate  at  high  frequencies. 

As  a  rule,  the  inductors  are  the  most  troublesome  elements.  The  quality 
of  resistors  and  capacitors  is  usually  much  higher  and  they  can  be  often 
treated  as  ideal,  particularly  at  moderate  frequencies. 

A  useful  equivalent  network  for  an  inductor  is  shown  in  Figure  45  [44] . 
This  three-element  network  offers  a  fairly  accurate  description  of  most  coils 
at  frequencies  where  they  are  designed  to  operate. 

As  an  example,  we  considered  the  effect  of  the  parasitic  elements  intro¬ 
duced  by  inductors  on  the  performance  of  the  equivalent  network  for  the 
dipole  antenna  (Figure  30) .  The  quality  factor  Q  =  200  at  frequency  w  =  10*10 
rad/sec  was  assumed  for  all  coils  and  the  computations  were  done  for  C  =  lpF 
and  C  =  2pF. 

The  network  responses  for  these  cases  are  compared  with  the  response  of 
the  original  network  in  Figures  46  and  47.  It  is  evident  that  the  parasitic 
effects  have  very  deteriorating  impact  on  the  network  performance.  As  can 
be  seen  from  these  figures  the  response  of  the  network  is  very  sensitive  to 
the  change  of  the  parasitic  capacitance.  Thus,  whenever  possible,  configura¬ 
tions  should  be  chosen  in  which  capacitors  appear  in  parallel  with  inductors. 
For  then  the  parasitic  effect  can  be  cancelled  by  reducing  the  size  of  the 
capacitor.  The  parasitic  capacitance  then  makes  up  for  the  reduction. 

5.3  NETWORK  SENSITIVITY 

Another  factor  limiting  the  practicality  of  the  realization  is  the  sensi¬ 
tivity  of  the  network  response  to  changes  in  element  values.  Thus,  it  is 


Equivalent  network  for  a  real 


Figure  46.  The  influence  of  the  parasitic  effects 


Figure  47.  The  influence  of  the  parasitic  effects: 
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important  to  choose,  from  all  the  networks  available,  the  one  which  exhibits 
the  least  sensitivity. 

It  was  shown  in  Section  TV  that  the  only  transformerless  synthesis  pro¬ 
cedure  available  for  minimum  biquadratic  functions  is  the  Bott-Duffin  method. 
Unfortunately,  the  Bott-Duffin  network  is  very  sensitive  to  component  errors. 
The  mechanism  behind  this  sensitivity  is  as  follows  [47]. 

The  admittance  of  the  Bott-Duffin  network  (Figure  21)  should  generally 
be  of  order  six.  Thus,  in  order  to  represent  a  second  order  function,  a 
cancellation  of  some  poles  and  zeros  must  occur  -  the  network  is  a  balanced 
bridge.  If  the  element's  values  are  changed  in  some  random  fashion,  then  all 
zeros  and  poles  migrate  and  the  bridge  balance  is  destroyed  producing  a 
significant  change  in  the  character  of  the  function. 

As  shown  in  Section  IV,  it  is  the  minimum  PR  character  of  the  approxi¬ 
mate  pole-pair  admittances  that  restricts  the  applicable  synthesis  tech¬ 
niques.  Thus,  it  is  conceivable,  that  by  excessive  resistive  padding,  but 
still  within  the  limits  of  acceptable  approximation,  other,  less  sensitive 
realization  could  be  obtained  (by  Miyata  procedure,  partial  fraction  expan¬ 
sion  ,  etc . ) . 

Finally,  it  should  be  noted  that  the  equivalent  networks  of  high-Q 
structures  (like  thin  wire  dipole  and  loop  antennas)  are  intrinsically  very 
sensitive  to  changes  in  element  values,  no  matter  what  the  form  of  their 
realization.  This  is  so  because  the  poles  of  the  sensitivity  function  coin¬ 
cide  with  the  poles  of  the  immittance  function. 


Section  VI 
CONCLUSIONS 

The  first  layer  conjugate  pole-pair  modules  for  the  dipole  and  loop  wire 
antennas  are,  in  general,  not  PR.  However,  they  can  be  made  PR  by  negligible 
resistive  padding.  It  is  conjectured  that  this  is  true  for  most  high-Q 
structures . 

The  padded  pole-pair  module  is  a  minimum-real-part  biquadratic  function. 
If  the  zero  of  the  real  part  of  this  function  occurs  at  zero  or  infinity, 
it  can  be  easily  realized  by  continuous  fraction  expansion  (Section  III) .  If 
however  the  zero  of  the  real  part  occurs  at  some  finite  frequency  (the  func¬ 
tion  is  Class  III,  in  terminology  of  Section  III)  then,  as  shown  in  Section 
IV,  the  Bott-Duffin  network  constitutes  the  exclusive  transformerless  form 
derivable. 

It  is  shown  that  for  structures  which  qualify  for  this  type  of  analysis 
the  Bott-Duffin  network  can  be  simplified  to  a  simple  four-element  network. 
This  network  has  been  shown  to  give  very  favorable  results  for  straight -wire 
and  wire- loop  antennas. 

Explicit  forms  of  this  and  other  networks  are  derived  in  Sections  III 
and  IV. 

It  is  shown  in  Section  V  that  the  early  time  response  of  the  dipole 
antenna  can  be  significantly  improved  by  capacitive  adjustment.  The  value 
of  the  capacitor  added  must  be,  however,  computed  apart  from  the  SEM 
representation . 

As  was  demonstrated  in  Section  V,  the  parasitic  elements  associated 
with  the  real  RLC  elements  can  have  a  very  deteriorating  influence  on  the 
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network  performance.  Therefore,  if  one  attempts  practical  realization  of  the 
network,  this  issue  should  be  carefully  considered  in  the  design. 


The  sensitivity  of  the  network  response  to  changes  in  element  values 
is  another  factor  which  limits  the  practicability  of  the  derived  networks. 
Thus,  whenever  possible,  the  least  sensitive  network  should  be  chosen  for 
realization. 

Finally,  the  equivalent  networks  of  high-Q  structures  are  intrinsically 
very  sensitive.  It  is  this  kind  of  sensitivity  which  dominates  in  the  case 
of  thin-wire  structures. 


L 


A 


Section  VII 


r 

* 

i 

is 


APPENDIX 

Below  the  listing  of  a  FORTRAN  program  for  the  conjugate  pole-pair  synthe¬ 
sis  of  an  admittance  function  implementing  the  results  of  Section  III  and 
Section  IV  is  included. 

The  program  works  in  conversational  mode.  The  input  consists  of  the 
poles,  the  residues,  and  the  normalization  constant  (identical  for  poles  and 
residues) .  The  program  computes  the  coefficients  of  the  biquadratic  (6) , 
the  minimum  real  part  frequency,  the  necessary  padding  (with  respect  to 
unmodified  pole-pair  module)  and  prints  out  the  element  values  (normalized 
and  absolute)  for  several  different  network  configurations.  The  element 
names  refer  to  relevant  figures  of  Sections  III  and  IV. 


A  sample  run  is  included  following  the  listing  of  the  program. 
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DIMENSION  ALFA (1 0) ,  BETA (10) ,  SIGMA  (10) , OMEGA  (1 0) 
COMMON  FAC 

987  CALL  INPUT(N, SIGMA, OMEGA,ALFA, BETA, FAC, SCALE) 

DO  1  I«1,N 
A=0. 

B=2.*ALFA(I) 

C-2 . * (ALFA ( I ) *S IGMA (I ) -BETA (I ) *OMEGA (I ) ) 

D=l. 

E=2. *SIGMA(I ) 

F-SIGMA (I ) **2+0 MEGA (I ) **2 
WRITE(5, 97)1 

97  FORMAT (///l OX' ******  POLE  PAIR  NO. '13'  ******'//) 

WRITE (5, 2 04) 

204  FORMAT ( / '  THE  COEFFICIENTS  BEFORE  PADDING') 
WRITE (5, 110)A,B,C,D,E,F 

110  FORMAT ( / '  A  ='E11.4'  B  -'Ell. 4'  C  -'Eli. 4/ 

1  '  D  -'Ell. 4'  E  -'Eli. 4'  F  -'Eli. 4/) 

CALL  FIND (A,B ,C ,D,E,F,W1,G1 ,IND ) 

Rl-O. 

IF(Gl.NE.O.  )R1— 1./G1 
WRITE (5, 205)W1 ,R1 

205  FORMAT (/'  THE  MINIMUM  REAL  PART  OCCURS  AT  THE 

1  (NORMALIZED)  '/'  FREQUENCY  -'E13.6/ 

1  '  THE  NECESSARY  PADDING  IS  -'E13.6'  OHMS'/) 
WRITE (5, 125) 

125  FORMAT ( /'  NEW  COEFFICIENTS  (AFTER  PADDING)') 

WRITE (5, 1 10)A,B,C,D,E»F 
IF (IND.EQ. 0)GOTO  51 
IF(IND.EQ.-1)CALL  TRY1 (A,B,C,D,E,F) 

IF(IND.EQ. 1)CALL  TRY2(A,B,C,D,E,F) 

GOTO  1 

51  CALL  BRUNE(A,B,C,D,E,F) 

CALL  HELP(A,B,C,D,E,F,KQ) 

IF(KQ.EQ. 0)CALL  ABOTT (A,B,C,D,E,F) 

IF(KQ.EQ. 1)CALL  BBOTT (A,B,C,D,E,F) 

I  CONTINUE 
WRITE(5,185) 

185  FORMAT (///'  ENTER  1  FOR  NEXT  PROBLEM 

1  OR  ZERO  TO  STOP'/) 

READ  (5,  *)KY 
IF(KY.EQ. 1)G0T0  987 
STOP 
END 
C 

SUBROUTINE  FIND(A,B,C,D,E,F,WO,GO,IND) 

LOGICAL  BUG 

V(X)-(W-X*U)/(2.*D*(A-X*D) ) 

BUG-. FALSE. 

ZO-C/F 

Zl-A/D 

IF(BUG)WRITE(5, 11)Z0,Z1 

II  FORMAT ( '  ZO  -'E13.5'  Z1  -'E13.5) 

Z-AMIN1 (Z0,Z1) 
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IND— 1 

Wl-0. 

IF (ZO.LE.Zl)GOTO  5 

IND-1 

W1-1.E38 

5  W«A*F+C*D-B*E 

U-2.*D*F-E*E 
A1 *U  *U-4. * (D*F ) **2 
B 1 -4 . *A*D  *F**2+4 . *C *D  **2  *F- 2 . *U  *W 
Cl«W**2-4. *A*C*D*F 
DEL«B1*B1-4.*A1 *C1 

IF(DEL.LT.O. )STOP'  FIND:  DELTA.LT.O  ?' 

SQ-SQRT(DEL) 

Xl*(-Bl-SQ)/2./Al 

X2-(-Bl+SQ)/2./Al 

Ol-V(Xl) 

02-V(X2) 

IF(BUG)WRITE(5,10)X1,X2,01,02 

10  FORMAT (  / '  XI  «'E13.5'  X2  ='E13.5/'  01  «'E13.5'  02  »'E13.5) 

IF(01.LT.O. .AND.02.LT. 0. )GOTO  7 
IF(Ol.LT.O. )X1“1.E38 
IF(O2.LT.0.)X2-l.E38 
IF(X1.LT.X2)G0T0  3 
G0=*X2 

WO=SQRT(02 ) 

GOTO  2 

3  G0=X1 

WO-SQRT(Ol) 

2  IF (GO .GE . Z )GOTO  7 

IND-0 

1  A-A-G0*D 

b-b-go*e 

C-C-G0*F 
RETURN 
7  GO-Z 

W0-W1 
GOTO  1 
END 
C 

SUBROUTINE  HELP(A,B,C,D, E,F,KQ) 

KQ-0 

AF-A*F 

Ql-SQRT(AF) 

CD-C*D 

Q2-SQRT(CD) 

BE-B*E 

Q3-SQRT(BE) 

Q4-Q1-Q2 

QQ-Q4**2-BE 

IF(ABS(QQ).GT.1.E-4)ST0P'  HELP:  NOT  A  MINIMUM  FUNCTION' 
IF(Q4.LT.0.)KQ«1 
RETURN 
END 
C 
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SUBROUTINE  ABOTT(A,B,C,D,E,F) 

REAL  M,L,L1,LN1,L2,LN2,LN 
COMMON  FAC 

M-B*SQRT  (D*F  )+E*SQRT  (A*C  ) 

CX-SQRT(A*B/E/F) 

L1-D/M*SQRT(D*F) 

Cl-M/D/F*SQRT(A/C) 

Rl-D/A 

L2-M/A/C*SQRT(A/C) 

K2-F/C 

C2-A*C/M*SQRT(D /F) 

L-D/C*SQRT(B*F/A/E) 

CN-CX/FAC 

LN1-L1/FAC 

CN1-C1/FAC 

LN2-L2/FAC 

CN2-C2/FAC 

LN-L /FAC 

WRITE (5, 131) 

131  FORMAT ( /'  BOTT-DUFFIN-CASE(A)-UNMODIFIED'/) 

WRITE (5, 132) 

132  FORMAT ( 7X'N0RMALIZED ' 6X'FINAL ' ) 

WRITE (5, 100)CX,CN,Ll ,LN1, Cl ,CNl ,Rl ,R1 ,R2, R2, 

1  L2,LN2,C2,CN2,L,LN 

100  FORMAT ( 

1  '  CO  -'2E13.5/'  LI  -'2E13.5/'  Cl  -'2E13.5/'  R1  -'2E13.5// 

2  '  R2  -'2E13.5/'  L2  -'2E13.5/'  C2  -'2E13.5/'  LO  -'2E13.5/) 
R-F*SQRT(B*D)4£*SQRT(C*E) 

CX«SQRT(A*B/E/F) 

L 1-R*D /M**2 *SQRT (B*D ) 

C 1 -M**2 /R/C /D  *SQRT (A*C /B /F ) 

Rl-D/A 

C2-A*C/R*SQRT(D*E/A/F) 

R2-F/C 

L-R/C/SQRT (C*E ) 

CN-CX/FAC 

LN1-L1/FAC 

CN1-C1/FAC 

CN2-C2/FAC 

LN-L/FAC 

WRITE(5, 133) 

133  FORMAT ( /'  BOTT-DUFFIN-CASE (A) -MODIFIED (A) '/) 

WRITE(5,132) 

WRITE (5 , 1 01 )CX, CN,L1 ,LN1 , Cl ,CN1,C2, CN2,L,LN,R1 , R1 ,R2 ,R2 

101  FORMAT C 

1  '  CO  -'2E13.5/'  LI  -'2E13.5/'  Cl  -'2E13.5/'  C2  -'2E13.5/ 

2  '  LO  -'2E13.5/'  R1  -'2E13.5/'  R2  -'2E13.5/) 

Q-A*SQRT (C*E)+B*SQRT (B*D ) 

CX-SQRT (A*B /E/F) 

CN-CX/FAC 

Rl-D/A 

C 1  -Q  /C  /D  *SQRT  ( A  *C  /B  /F ) 

CN1-C1/FAC 

L2-M**2/Q/C/SQRT(E*C) 
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LN2-L2/FAC 

L»D/Q*SQRT(B*D) 

LN-L/FAC 
R2-F/C 
WRITE(5 ,134) 

134  FORMAT (/'  BOTT-DUFFIN-CASE (A) -MODIFIED (B) '/) 

WRITE(5, 132) 

WRITE (5, 102)CX,CN,C1,CN1,L2,LN2,C2,CN2,L,LN,R1 ,R1,R2,R2 
102  FORMAT ( 

1  '  CO  -'2E13.5/'  Cl  -'2E13.5/'  L2  -'2E13.5/'  C2  -'2E13.5/ 

2  '  LO  -'2E13.5/'  R1  -'2E13.5/'  R2  -'2E13.5/) 

CX-S  QRT (A*B / E /F ) 

UN-CX/FAC 

Rl-D/A 

C 1-Q**2*SQRT (A/C ) / (M*B*D ) 

CNl-Cl/FAC 

L1-M*B*D/ (C*Q*SQRT (A*E) ) 

LN1-L1/FAC 

L-D/Q*SQRT(B*D) 

LN-L/FAC 

L2-M*SQRT (A/C ) / <A*C ) 

LN2-L2/FAC 

C2-A*C*SQRT(D/F)/M 

CN2-C2/FAC 

R2-F/C 

WRITE (5, 152) 

152  FORMAT (/'  DARLINGTON-CASE(A) '/) 

WRITE (5, 132) 

WRITE(5,100)CX,CN,L1,LN1,C1,CN1,R1,R1,R2,R2, 

1  L2,LN2,C2, CN2,L ,LN 

RETURN 
END 
C 

SUBROUTINE  BRUNE(A,B,C,D,E,F) 

REAL  LA,LB,M,LNA,LNB,MN 

COMMON  FAC 

LA-D/B 

LNA-LA/FAC 

LB-A*F/B/C 

LNB-LB/FAC 

M-A/B*SQRT(D*F/A/C) 

MN-M/FAC 

CX-B/F 

CN-CX/FAC 

R-F/C 

WRITE(5, 135) 

135  FORMAT (/'  BRUNE  NETWORK'/) 

WRITE(5,132) 

132  FORMAT (7X'N0RMALIZED ' 6X'FINAL ' ) 

WRITE (5  > 105)LA,LNA,LB,LNB ,M,MN , CX, CN,R,R 
105  FORMAT ( 

1  '  LA  -'2E13.5/'  LB  -'2E13.5/'  M  -'2E13.5/ 

2  '  CO  -'2E13.5/'  R  -'2E13.5/) 

RETURN 
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END 
C 

SUBROUTINE  TRY1 (A,B,C,D,E,F) 

REAL  Ll.LNl 
COMMON  FAC 

IF (ABS (C) «GT. l.E-8) STOP'  TRY  1:  C.NE.O  ?' 

Q-B*E-A*F 

IF(Q.LT.O. )STOP'  TRY  1:  BE-AF  .LT.  ZERO  ?' 

IF(Q/B**2.GT.D/A)G0T0  7 

Cl-B/F 

CN1-C1/FAC 

R1»Q/B**2 

Ll«D/B-A*Rl/B 

LN1-L1/FAC 

R2-B*L1/A 

WRITE(5, 107) 

107  FORMAT ( /  '  LADDER  NETWORK  1-A  (ZERO  AT  ZERO)'/) 

WRITE (5, 132) 

132  FORMAT (7X'NORMALIZED'6X' FINAL') 

WRITE (5 , 120)C1,CN1,R1,R1,L1,LN1,R2,R2 

120  FORMAT ( '  Cl  -'2E13.5/'  R1  -'2E13.5/'  LI  -'2E13.5/ 

1  '  R2  «'2E13.5/) 

RETURN 
7  Cl-B/F 

CN1-C1/FAC 

R1*=D/A 

C2«A**2*B/ (A*Q-B**2*D) 

CN2-C2/FAC 
R2-Q/B**2-D/A 
WRITE(5 ,137) 

137  FORMAT ( / '  LADDER  NETWORK  1-B  (ZERO  AT  ZERO)'/) 

WRITE (5, 132) 

WRITE (5, 121)Cl,CNl,Rl,Rl , C2,CN2,R2,R2 

121  FORMAT ( '  Cl  -'2E13.5/'  R1  -'2E13.5/'  C2  -'2E13.5 / 

1  '  R2  -'2E13.5/) 

RETURN 

END 

C 

SUBROUTINE  INPUT(N, SIGMA, OMEGA, ALFA, BET A, FAC, SCALE) 
REAL  ALFA(IO) ,BETA(10) ,SIGMA(10) ,OMEGA(10) 

WRITE(5, 200) 

200  FORMAT (//'  GIVE  N  -  THE  NUMBER  OF  CONJUGATE 
1  POLE  PAIRS'/) 

READ (5, *)N 
WRITE(5,201) 

201  FORMAT ('  GIVE  THE  POLES  AND  THE  RESIDUES'/ 

1'  -REAL  +  j  IMAG  4REAL  +  j  IMA G'/) 

DO  37  I-l.N 

37  READ (5, *) SIGMA (I) , OMEGA (I ) , ALFA (I) ,BETA(I ) 

WRITE(5,203) 

203  FORMAT (/'  GIVE  THE  NORMALIZATION'/) 

READ (5, *)FAC  . 

RETURN 

END 
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C 

SUBROUTINE  TRY2 (A,B,C,D,E, F) 

REAL  L1,LN1,L2,LN2 
COMMON  FAC 

IF (ABS (A) .GT. 1. E-8)ST0P'  TRY  2:  A.NE.ZERO  ?' 
q«B*E-C*D 

IF(Q.LT.O. )STOP'  TRY  2:  BE-CD  .LT.  ZERO  ?' 

IF(Q/B**2.GT.F/C)G0T0  17 

Ll-D/B 

LN1-L1/FAC 

Rl-Q/B**2 

C1«B**3/(B**2*F-Q*C) 

CN1-C1/FAC 
R2-F/C-Q/B**2 
WRITE (5, 108) 

108  FORMAT (/'  LADDER  NETWORK  2-A  (ZERO  AT  INFINITY)'/) 

WRITE (5, 132) 

132  FORMAT(7X'NORMALIZED'6X'FINAL') 

WRITE (5,1 20 )C1,CN1,R1,R1 ,L1,LN1,R2,R2 
120  FORMAT ( '  Cl  -'2E13.5/'  R1  -'2E13.5/'  LI  -'2E13.5/ 

1  '  R2  «'2E13.5/) 

RETURN 
17  Ll-D/B 

LN1-L1/FAC 

Rl-F/C 

L2«(Q*C-B**2*F)/(B*C**2) 

LN2-L2/FAC 
R2-Q/B**2-F/C 
WRITE(5, 138) 

138  FORMAT (/'  LADDER  NETWORK  2-B  (ZERO  AT  INFINITY)'/) 

WRITE (5, 132) 

WRITE(5, 123)L1,LN1,R1 ,Rl ,L2,LN2,R2,R2 
123  FORMAT ( '  LI  -'2E13.5/'  R1  -'2E13.5/'  L2  -'2E13.5/ 

1  '  R2  -'2E13.5/) 

RETURN 

END 

C 

SUBROUTINE  BBOTT(A,B,C,D,E,F) 

REAL  M,L,LN,L1,LN1,L2,LN2 
COMMON  FAC 

M-B*SQRT(D*F)+E*SQRT(A*C) 

CX-A/F*SQRT(C*E/B/D) 

CN-CX/FAC 
L1-M/A/C*SQRT (D/F) 

LN1-L1/FAC 

C1-A/M*SQRT(A*C) 

CN1-C1/FAC 

Rl-D/A 

C2«M/F/SQRT(D*F) 

CN2-C2/FAC 

L2-D*F/M*SQRT(A/C) 

LN2-L2/FAC 

R2-F/C 

L-SQRT(D *E/B/C) 
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LN-L/FAC 
WRITE(5, 231) 

231  FORMAT (/'  BOTT-DUFFIN  CASE(B ) -UNMODIFIED' /) 
WRITE (5, 132) 

132  FORMAT(7X'NORMALIZED'6X'FINAL' ) 

WRITE (5, 200)CX, CN,Ll»LNl ,C1,CN1,R1,R1, 

1  C2,CN2,L2,LN2,R2,R2,L,LN 

200  FORMAT ( 

1  '  CO  -'2E13.5/'  LI  -'2E13.5/'  Cl  -'2E13.5/'  R1  - 

2  '  C2  -'2E13.5/'  L2  -'2E13.5/'  R2  -'2E13.5/'  LO  - 
Q-D*SQRT(B*F)+E*SQRT(A*E) 

CX-A/Q*SQRT (A*E) 

CN-CX/FAC 

Rl-D/A 

L1-Q/A/C*SQRT(C*D/E/F) 

LNl-Ll/FAC 

R2-F/C 

C2-M**2/Q/F/SQRT (B*F ) 

CN2-C2/FAC 

L2=Q*A*F/M**2*SQRT(B*D/A/C ) 

LN2-L2/FAC 

L«SQRT(D*E/B/C) 

LN-L/FAC 
WRITE(5, 233) 

233  FORMAT ( / '  BOTT-DUFFIN-CASE (B )-MODIFIED (A) ' / ) 
WRITE (5, 132) 

WRITE (5,201)CX,CN,R1,R1 ,L1,LN1 ,R2,R2, 

1  C2,CN2,L2,LN2,L,LN 

201  FORMAT  ( 

1  '  CO  -'2E13.5/'  Rl  -'2E13.5/'  Ll  -'2E13.5/'  R2 

2  '  C2  -'2E13.5/'  L2  -'2E13.5/'  LO  -'2E13.5/) 
R-C*SQRT(A*E)«*SQRT(B*F) 

CX-R/F/SQRT(B*F) 

CN-CX/FAC 

L2-A*F/R*SQRT(B*D/A/C) 

LN2-L2/FAC 

L1-M**2/R/A/C*SQRT(C*D/E/F) 

LNl-L 1/FAC 

C1-R*A/M**2*SQRT(A*E) 

CNl-Cl/FAC 

Rl-D/A 

R2-F/C 

L-SQRT (D  *E /B /C ) 

LN-L/FAC 
WRITE(5  234) 

234  FORMAT (/'  BOTT-DDFFIN-CASE (B) -MODIFIED (B) '/) 
WRITE(5,132) 

WRITE(5,202)CX,CN,L2,LN2,L1,LN1,C1,CN1, 

1  R1,R1,R2,R2,L,LN 

202  FORMAT ( 

1  '  CO  -'2E13.5/'  L2  -'2E13.5/'  Ll  -'2E13.5/'  Cl 

2  '  Rl  -'2E13.5/'  R2  -'2E13.5/'  LO  -'2E13.5/) 

CX-R /F /SQRT (B*F ) 

CN-CX/FAC 


'2E13.5// 

'2E13.5/) 


'2E13.5/ 


-'2E13.5/ 


L2-M*B*F*SQRT(A/C ) /R**2 
LN2-L2/FAC 
L1«M*SQRT(D/F)/A/C 
LN1-L1/FAC 

C2»R*A*SQRT(C*E)  /M/B/F 

CN2-C1/FAC 

Rl-D/A 

R2-F/C 

L«SQRT(D*E/B/C) 

LN-L/FAC 

C1-A*SQRT(A*C)/M 

CN1-C1/FAC 

WRITE(5,252) 

FORMAT (/'  DARLINGTON-CASE(BOV) 

WRITE (5, 132) 

WRITE(5,200)CX,CN,L1,LN1,C1,CN1,R1,R1, 
1  C2,CN2,L2,LN2,R2,R2,L,LN 

RETURN 
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ex  granit.for 
LINK:  Loading 

[LNKXCT  GRANIT  Execution] 


GIVE  N  -  THE  NUMBER  OF  CONJUGATE  POLE  PAIRS 


2 


GIVE  THE  POLES  AND  THE  RESIDUES 
-REAL  +  j  IMAG  +REAL  +  j  IMAG 

.08427  .9158  .1112e-2  .3121e-3 
.1473  2.870  . 1319e-2  .3301e-3 

GIVE  THE  NORMALIZATION 

9.4248e6 


******  POLE  PAIR  NO.  1  ****** 


THE  COEFFICIENTS  BEFORE  PADDING 

A  -  O.OOOOE+OO  B  -  0. 2224E-02  C  — 0.3842E-03 
D  -  0.1000E+01  E  -  0.1685E+00  F  -  0.8458E+00 


THE  MINIMUM  REAL  PART  OCCURS  AT  THE  (NORMALIZED) 
FREQUENCY  -  0.000000E+00 
THE  NECESSARY  PADDING  IS  -  0.220129E+04  OHMS 


NEW  COEFFICIENTS  (AFTER  PADDING) 

A  -  0.4543E-03  B  -  0.2301E-02  C  -  0.0000E+00 
D  -  0. 1000E+01  E  -  0.1685E+00  F  -  0.8458E+00 


LADDER  NETWORK  1-A  (ZERO  AT  ZERO) 


NORMALIZED 
Cl  -  0.27200E-02 

R1  -  0.66342E+00 

LI  -  0. 43454E+03 
R2  -  0.22006E+04 


FINAL 

0.28860E-09 

0.66342E+00 

0.46107E-04 

0.22006E+04 
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******  pole  PAIR  NO.  2  ****** 


THE  COEFFICIENTS  BEFORE  PADDING 

A  -  0.0000E+00  B  -  0.2638E-02  C  — 0.1506E-02 
D  -  0.1000E+01  E  -  0.2946E+00  F  =*  0.8259E+01 


THE  MINIMUM  REAL  PART  OCCURS  AT  THE  (NORMALIZED) 
FREQUENCY  -  0.I59812E+01 
THE  NECESSARY  PADDING  IS  -  0.495870E+04  OHMS 


NEW  COEFFICIENTS  (AFTER  PADDING) 

A  -  0. 2017E-03  B  -  0.2697E-02  C  -  0.1593E-03 
D  «  0. 1Q00E+01  E  -  0. 2946E+00  F  -  0.8259E+01 


BRUNE  NETWORK 


NORMALIZED 

FINAL 

LA  - 

0.37073E+03 

0.39335E-04 

LB  - 

0.38764E+04 

0.41130E-03 

M  - 

0. 11988E+04 

0.12720E-03 

CO  - 

0.32662E-03 

0.34655E-10 

R  - 

0. 5I850E+05 

0.51850E+05 

BOTT-DUFFIN-CASE (A) -UNMODIFIED 


NORMALIZED 

FINAL 

CO 

m 

0.47285E-03 

0.50170E-10 

LI 

m 

0.36822E+03 

0.39069E-04 

Cl 

m 

0. 10634E-02 

0. 1 1283E-09 

R1 

m 

0.49587E+04 

0.49587E+04 

R2 

m 

0.51850E+05 

0.51850E+05 

L2 

m 

0.27340E+06 

0.29008E-01 

C2 

m 

0, 14322E-05 

0. 15196E-12 

LO 

m 

0.12157E+06 

0.12899E-01 

BOTT-DUFFI N-C ASE ( A ) -MOD IF I ED ( A ) 


NORMALIZED 

FINAL 

CO  - 

0.47285E-03 

0.50170E-10 

LI  - 

0.36745E+03 

0. 38987E-04 

Cl  - 

0. 10656E-02 

0. 1 1306E-09 

C2  * 

0. 99133E-06 

0.10518E-12 

LO  - 

0.39497E+06 

0.41908E-01 

R1  * 

0.49587E+04 

0.49587E+04 

R2  - 

0.51850E+05 

0.5I850E+05 

BOTT-DUFFIN-CASE (A )-MODIFIED (B ) 

NORMALIZED 

FINAL 

CO  = 

0.47285E-03 

0.50170E-10 

Cl  - 

0. 10666E-02 

0.1 1317E-09 

L2  - 

0.39460E+06 

0.41869E-01 

C2  - 

0.99133E-06 

0. 10518E-12 

LO  - 

0.36711E+03 

0.38951E-04 

R1  - 

0.49587E+04 

0.49587E+04 

R2  - 

0.51850E+05 

0.51850E+05 

DARLINGTON-CASE (A) 


NORMALIZED 

FINAL 

CO  ** 

0.47285E-03 

0.501 70E-10 

LI  * 

0.12121E+06 

0.12860E-01 

Cl  - 

0.10698E-02 

0.11351E-09 

R1  - 

0.49587E+04 

0.49587E+O4 

R2  - 

0.51850E+05 

0.51850E+05 

L2  - 

0.27340E+06 

0.29008E-01 

C2  = 

0. 14322E-05 

0.15196E-12 

LO  - 

0.36711E+03 

0.38951E-04 

ENTER  1  FOR  NEXT  PROBLEM  OR  ZERO  TO  STOP 

0 

STOP 

END  OF  EXECUTION 

CPU  TIME :  0.38  ELAPSED  TIME:  3:15.70 
EXIT 
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